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c/o  Dr.  Thomas  Kinder 
Office  of  Naval  Research 
800  N  Quincy  Street 
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Dear  Sir/Madam: 

This  letter  constitutes  the  revised  final  report  for  ONR  contract  N00014-88-K-0203 
awarded  to  Old  Dominion  University  for  the  period  October  i,  1987  though  January 
31, 1989.  The  original  final  report  was  submitted  to  the  Office  of  Naval  Research  as 
part  of  a  renewal  proposal. 

The  bulk  of  the  research  activities  conducted  under  this  contract  consisted  of  basic 
research  on  ocean  flow  dynamics  as  it  pertains  to  the  prediction  of  ocean  motion. 
During  the  contract  period,  three  papers  describing  aspects  of  the  research  were 
published  in  the  adjudicated  scientific  literature.  These  papers  are: 

1.  “Genesis  of  the  Gulf  of  Mexico  Ring  as  Determined  from  Kinematic  Analysis,” 
J.  Geophys.  Res.,  92(Cll),  11727-11740,  1987. 

2.  “Observed  and  Simulated  Kinematic  Properties  of  Loop  Current  Rings,”  J.  Geo¬ 
phys.  Res.,  93(C2),  1189-1198,  1988. 

3.  “Notes  on  the  Cluster  Method  for  Interpreting  Relative  Motions,”  J.  Geo¬ 
phys.  Res.,  93(C8),  ^  )7-9339,  1988. 

In  addition  to  these  papois,  the  results  of  the  research  were  reported  at  the  1988  Liege 
Colloquium  on  Ocean  Hydrodynamics,  Mesoscale/Synoptic  Coherence  in  Geophysical 
Turbulence.. 

As  a  part  of  the  research  effort,  the  funds  were  used  to  support  William  Indest,  a 
Ph.D.  graduate  student  in  the  Department  of  Oceanography. 

I  am  grateful  to  the  Office  of  Naval  Research  for  continued  support  for  this  research. 
Respectfully, 

A.  D.  Kirwan,  Jr. 

Samuel  L.  and  Fay  M.  Slover 
Chair  of  Oceanography 
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Genesis  o&a  Gulf  of  Mexico  Ring  as  Determined  From  Kinematic  Analyses 

Jamcs  K.  Lewis 

Srkpee  .IfpUe^inM  Immminmid  Ctrfermiem.  CeKrqt  Stmim,  Texts 

D.  Kirwan.  Jr.' 

Uturmilr  uf  Sitth  TUnla.  Orftrtmrm  tf  ilarime  Sdesee.  St.  TeterAmg 

The  kiiKiiutics  of  ihe  Loop  Currcni  m  tio^ei  usmg  tniectories  of  dtihm  to  the  Guir  of  Mrtim 
during  mid-June  through  Septemiief  I9i5.  One  of  the  driften  was  in  the  Loop  Cumat  proper,  while 
iilhcr  driften  were  in  Inm  reocnily  shed  tjintp  Current  tinps.  The  drifter  in  the  Loop  Oirtent  ihnwed 
5ir<m§  aniicychmic  mulion  during  the  fludy  period.  Uiir  Loop  Cumenl  anticyclone  first  began  oK  the 
northwest  coast  of  Cuba.  It  rapidly  moned  nerthwrard  into  the  Guif  of  Meaico  as  a  ring  pinched  ofiTftoni 
the  Loop  Current.  Analysis  of  the  Loop  Current  drifter  motion  showed  that  the  anlicydone  became  an 
integral  part  of  the  Loop  Currents  taking  on  many  of  the  characteristics  of  the  most  recently  shed  ring. 

The  residts  of  the  .sn.slysis  suggest  a  process  by  which  Loop  Current  rings  can  be  generated.  Apparently, 
this  mechanism  on  cause  the  Loop  Current  to  become  reconfigured  in  2-3  nionihs  for  beginning  the 

process  of  ring  scpariilHin. 


I.  lNrR(M)t.'(TK>N 

ITie  shedding  of  Loop  Current  rings  has  a  major  impact  on 
processes  in  the  central  and  western  Gulf  of  Mexico.  These 
large  aniicyelones  transport  a  tremendous  amount  of  momen¬ 
tum.  heat,  and  salt  across  the  gulf,  all  the  way  to  the  Mexican 
coast  [Elliot.  1982:  Kirwan  el  al.,  1984;  Lewis  and  Kirwan, 
19X5  j.  In  <»r<lcr  t<>  consider  halanccs  of  momenttim.  muvs.  and 
heat  within  the  Gulf  of  Mexico,  it  is  important  to  have  some 
icic.i  <if  the  characteristics  of  the  kinematics  of  the  Loop  Cur¬ 
rent.  including  how  often  a  ring  may  be  pinched  olT. 

During  June  1985.  an  attempt  was  made  to  put  an  Argos 
drifter  into  the  Loop  Current  as  a  ring  was  pinching  oflT.  How¬ 
ever,  the  ring  did  not  totally  disconnect  from  the  l.oop  Cur¬ 
rent.  and  the  drifter  ended  up  spending  approximately  3 
months  in  the  Loop  Current  proper.  The  drifter  exited  the 
Gulf  of  Mexico  (GOM)  through  the  Florida  Straits  in  Septem¬ 
ber  1985  (Figure  I.  drifter  ,t.354).  Fortunately,  a  second  drifter 
(.3378.  Figure  I)  was  placed  in  the  ring  in  July  1985.  Together, 
these  two  Lagrangian  data  sets  provide  a  unique  view  of  the 
motion  uf  the  Loop  Current  while  it  is  in  the  process  of  shed¬ 
ding  a  ring.  They  provide  special  insight  into  Loop  Current 
kinematics  as  this  system  extends  northward,  as  well  as  the 
kinematics  of  a  Loop  Current  ring.  Finally,  we  comment  on 
the  existence  of  a  third  anticyclone  in  the  GOM.  The  eddy 
was  di.scuvered  sererdipitously  as  a  third  drifter  became  en¬ 
trained  in  its  (low  held  and  circled  around  the  eddy  all  the 
way  to  the  Mexican  coast  (G.  Forrisiall.  personal  communi¬ 
cation.  1986).  Thi.v  eddy  Inis  been  referred  to  as  Ghost  Eddy 
because  (here  had  been  no  previous  indication  of  its  existence 
prior  to  the  first  of  August  1985.  Its  general  motion  and  cxicnl 
during  August  and  September  1985  is  shown  in  Figure  I. 
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The  following  analysis  is  obtained  from  the  movement  of 
drifters  3354  and  3378  as  well  as  concurrent  sea  surface  tm- 
perature  (SST)  data  and  XBT  data.  These,  along  with  the 
presence  and  location  of  Ghost  Eddy,  suggest  a  new  process 
by  which  an  anlicyclonic  vortex  is  formed  that  can  eventually 
become  a  Loop  Current  ring.  The  original  mechanism  for 
spinning  up  the  vortex  appears  to  be  the  lateral  sharing 
stresses  caused  by  the  Loop  Current  olT  the  northwest  coast  of 
Cuba. 

A  detailed  description  of  the  kinematic  characteristics  of  the 
Loop  Current  and  the  ring  is  provided  by  an  analysis  of  the 
drifter  paths.  This  analysb  gives  us  the  time  histories  of  the 
rotation  r.iie,  cccentria'ty  and  orientation  of  the  ellipses  of  the 
tr.-rjectories.  swirl  -velocities,  and  movement  of  the  centers  of 
rotation.  A  comparison  is  made  of  the  kinematics  of  the  Loop 
Current  as  delemiined  by  drifter  3354  during  the  time  it  co¬ 
existed  with  a  ring.  The  latter  is  designated  by  drifter  3378. 
The  comparison  of  the  last  part  of  drifter  3354  (before  it  exited 
the  GOM)  with  the  first  month  of  ring  3378  provides  an  in¬ 
teresting  contrast  between  the  motion  characteristics  of  a 
northern  extension  of  the  Loop  Current  before  and  after  it 
pinches  off  to  become  a  ring. 

2.  Path  Data 

The  path  data  used  in  this  study  are  the  position  data  of  the 
drifters  with  Argos  identification  numbers  3354  and  3378.  The 
GOM  ring  that  was  seeded  by  drifter  3378  will  be  referred  to 
as  ring  3378.  References  to  drifter  3354  will  indicate  a  refer¬ 
ence  to  the  Loop  Current  proper.  Drifter  3378  was  drogued  by 
a  weighted  200-m  line,  while  drifter  3354  was  drogued  by  a 
window  shade  drogue  at  100  m. 

Drifter  3354  was  seeded  at  25.9“N.  87.9°W  on  June  18, 1985. 
The  drifter  immediately  moved  .southeastward  ~525  km  and 
reached  23.2’'N.  83.7'W  by  June  29,  1985.  At  this  point  the 
drifter  became  entrained  in  a  westward,  aniicyclonic  flow  field 
with  a  center  of  rotation  al  about  24°N,  85.5°W.  After  two 
rotations  (into  mid-August)  the  drifter  suddenly  moved  north¬ 
westward  and  made  three  additional  rotations  centered  at 
about  25,5'  N,  Sfi.S’W.  The  drifter  then  left  this  flow  pattern  in 
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Fig.  I.  Trajectories  of  drifter  3.^54  from  tnid-June  through  mid-September  ITUS  and  of  drifter  3378  from  mid-July 
through  September  1985.  Also  shown  is  the  lociiion  of  Ghost  Eddy  duiing  August  and  September  I98S.  Squares  denote 
the  beginning  positions  of  the  drifter  trajectories  and  triangles  denote  tlic  end  positions. 


niid-Scpiember  1985  and  exiled  the  GOM  through  the  Flori¬ 
da  Straits. 

Drifter  3.378  was  .seeded  in  ring  .3378  at  26.4  -N.  H9..3  W  on 
July  IX.  198.5.  Ily  that  lime,  the  ring  had  totally  separated 
front  the  l.oop  t'urreni.  As  shown  in  l'■ipure  I.  ring  3378 
.slowly  mo\ed  westward,  reaching  approxiinaiely  the  91  W 
meridian  by  mid-August.  Recall  that  after  mid-August,  the 
Loop  Current  (as  determined  by  drifter  ,3.3.54)  mosed  north- 


TABI.n  I.  Listing  of  Cfiiises  During  Which  .XBT  Data  Were 
Collected  in  the  C'cntr.il  and  Easiern  Gulf  of  Mexico 


Vessel 

DiiietlWi 

F.  M.  Qmriiv 

May  7  X 

M,  V  1 

Mav  to  1 1 

M  V  SicttU  Uhpmiiu 

May  17  18 

K.  St.  Qiicrnr 

May  26  27 

M/V  Sii’im  //(spimm 

May  27  28 

M'V  Ithpama 

May  .«K3I 

E.  SI.  Qiiceiiy 

June  7  8 

L  M.  Qmviiy 

June  I.3-I4 

M,'V  Swim  tli^puniii 

June  26-28 

M/V  Swim  IlixpiiiiM 

June  29-.30 

E.  St.  Qiu'i'ny 

July  1-2 

R/V  SiiiiriMswr 

July  9  22 

M,'V  Not  Co  6 

July  16-19 

MW  Sii'iui  lli'ipiiiuii 

July  16-17 

E.  SI.  Qiiirnv 

July  31  to  Aug.  1 

MW  Swim  llispoiihi 

Aug.  16-18 

MW  Aiiibiixxiiilor 

Sept.  4  S 

M,  V  Aiiihiixmliir 

Sept.  13  14 

M/V  AiiihiisuiJiir 

Sept.  24-2.5 

westward  and  remained  there  until  at  least  mid-September. 
During  that  same  JOnfay  period,  ring  .3378  continued  moving 
westward,  reaching  to  the  92.5'W  meridian. 

The  onshore  industry  launched  a  drifter  in  June  1985  which 
eventually  became  entrained  in  the  flow  held  of  Ghost  Eddy. 
At  the  beginning  of  Atigu.si  1985.  this  ring  was  centered  at 
23.5  N.  93  W  (Figure  I).  The  ring  gradually  moved  southwest- 
ward  across  the  deepest  portion  of  the  GOM.  By  the  end  of 
September  1985.  Ghost  Eddy  was  still  strongly  rotating  at 
23  N.  94.5  W. 


3.  TkMW-RATURK  DaTA 

A  number  of  XBT  data  .sets  were  collected  in  the  eastern 
GOM  during  May-Seplember  1985  (Table  1).  These  data  pro¬ 
vide  indications  of  the  location  and  the  structure  of  the  Loop 
Current  and  ring  3378.  In  addition,  weekly  SST  contour 
charts  arc  available  for  the  entire  gulf.  Some  of  these  SST  data 
are  presented  here,  but  unfortunately  little  in  the  way  of  ther¬ 
mal  structure  cun  be  picked  out  .since  the  SST  gradients  are 
small  in  the  GOM  during  the  study  period  (summer  and  early 
fall). 

Cruise  tracks  for  the  period  May  26-31,  1985.  arc  shown  in 
Figure  2ii.  The  resulting  XBT  data  are  shown  in  Figures 
2/)  2tt.  These  data  were  collected  before  ring  3378  was  shed, 
and  the  edges  of  the  Loop  Current  (maximum  horizontal  tem¬ 
perature  gradients)  indicated  by  (he  XBT  data  are  shown  in 
Figure  2ii.  (Smaller  slopes  of  the  isotherms  were  taken  as  an 
indication  that  the  transects  were  crossing  the  Loop  Current 
edge  at  tin  tingle,  and  some  of  the  edges  were  thus  drawn  at  a 
45‘  angle  to  the  cruise  track.)  It  appears  that  the  Loop  Cur- 
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Fig.  2.  (a)  Cruiw  irscks  Tor  XBT  dais  collected  during  May  26'3I.  I9R5.  The  arrows  denote  the  flow  at  the  edges  of 
the  Loop  Current  based  on  the  vertical  temperature  structure  shown  in  {h)  temperature  data  from  the  £.  M.  Quttny  cruise, 
May  2t>-27, 198.\  (cl  lempcraiure  data  from  the  SifiM  llhpania  cruise.  May  27-211, 1985,  and  (d|  temperature  data  from  the 
SieiM  Hispaiila  cruise.  May  30'31. 1985. 


rent  w.is  flowing  northward  on  the  east  side  of  the  Yucatan 
Straits,  swirled  westward  and  around  the  region  where  ring 
.3378  would  be  spawned,  and  then  turned  northeastward 
before  exiting  the  COM. 

As  previously  mentioned,  drifter  3354  was  deployed  in  the 
northward  extension  of  the  Loop  Current  during  mid-June 
I985.  The  trajectory  of  3354  during  June  I8-24.  I985.  is  plot¬ 
ted  on  the  corresponding  COM  SST  chart  in  Figure  3.  As 


indicated  by  the  trajectory  (and  to  some  degree  by  the  SST 
data),  ring  3378  had  obviously  not  separated  from  the  Loop 
Current  during  this  period.  The  return  flow  from  the  area  of 
the  northward  extension  appears  to  be  mostly  southeastward, 
toward  the  northern  coast  of  Cuba. 

Three  XBT  data  sets  were  collected  during  June  26  to  July 
2.  I985  (Figure  4),  and  the  temperature  profiles  are  shown  in 
Figures  4h-4J.  The  outline  of  ring  3378  is  quite  distinct  in 
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Fig.  3.  Tfjijtclory  of  drifter  3J54  and  SST  data  for  June  18-24. 1985. 


I'igurcs  4/i  and  4f,  while  the  Loop  Current  can  be  .seen  to 
estend  to  25  N  in  Figure  4ii.  Tliesc  data  would  seem  to  imply 
that  ring  .3.t7X  is  still  connected  to  the  Loop  Current.  The  SST 
data  (or  June  29  to  July  5.  1985.  along  with  the  associated 
trajectory  of  drifter  3.354  arc  shown  in  Figure  5.  The  outlines 


of  the  ring  and  the  Loop  Current  from  the  XBT  data  in 
Figure  4<i  arc  also  shown,  and  we  see  some  of  the  rotational 
characteristics  within  the  Loop  Current  from  the  path  of  drif¬ 
ter  .3354.  The  rotation  near  the  Cuban  coast  continued  into 
mid-July. 
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Fig.  4.  (ii)  Cruise  tracks  for  Xltt'  dalii  cullecicd  during  June  26  tu  July  2.  1985.  Uie  arrows  denote  the  flow  at  the 
edges  of  the  t.oop  Current  hased  on  the  vertical  Icmpcriittirc  structure  shown  in  (h)  temperature  data  from  (he  Sieiui 
IlhponUi  cruise.  June  26-28.  1985.  (c)  temperature  data  from  the  Siena  llhpania  cruise,  June  29-30,  1985,  and  (d) 
(eniperature  data  from  (he  K.  M.  Qmviiy  cruise.  July  1-2. 1985. 
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During  July  16  19.  1985.  (wo  XBT  Jala  sets  were  collected 
(Figure  6)  which  distinctly  show  that  ring  3378  had  separated 
froin  the  Loop  Current.  The  S.ST  contours  for  the  .same  period 
are  shown  in  Figure  7.  There  is  no  surface  signature  of  ring 
3378.  although  its  rotation  is  well  delineated  by  the  trajectory 
of  drifter  3378.  As  for  the  Loop  Current.  Figure  7  shows  that 
the  rotational  feature  within  (he  current  had  reached  the 
northern  edge  of  the  current,  about  25‘N. 

Both  (he  l.cKip  Current  and  ring  3378  continued  to  rotate 
anticrclonicaliy  through  August  -1985.  An  XBT  data  .set  col¬ 
lected  on  August  16-18.  1985  (Figure  8)  indicates  that  the 
Lot'p  Current  had  e.xtendcd  northward  to  26.6  N.  Looking  at 
the  S.ST  map  for  August  12-19.  1985  (Figure  9).  wc  see  that 
ring  L^7S  and  the  Loop  Current  were  rotating  at  ahrmt  (he 
same  latitude,  25.5  N.  By  inid-Septenihcr,  (he  sea  .surface  had 
cooled  suflicicntly  in  the  eastern  (JOM  .so  that  the  northern 


extension  of  the  Loop  Current  is  readily  defined  (Figure  10). 
The  trajectories  for  this  period  of  time  show  ring  3378 
strongly  rotating  but  drifter  3354  leaving  the  rotational  fea-. 
ture  of  the  Loop  Current  along  the  eastern  side  of.  a  28°C 
tongue  of  water.  The  following  week's  SST  chart  (Figure  II) 
shows  the  28  C  tongue  reaching  up  to  27.5”N,  while  drifter 
33.54  moved  toward  the  Florida  Straits  along  the  northern 
coast  of  Cuba.  Ring  3378  was  still  rotating  with  a  center  at 
about  25.5  N.92  W. 

4.  Kinoiatic  Analyses 

The  trajectory  data  were  used  to  calculate  various  kin¬ 
ematic  parameters  of  the  flow  ticid.  The  methodology  outlined 
by  Kirwan  tt  al.  [1984.  1987]  for  geophysical  fluids  was  used 
in  thc.se  calculations,  (l-’or  more  details,  see  Kirwan  tt  al. 
[1987].)  Aillunigh  Kirwon  tt  o/.  [1984]  used  a  lOO-hqur  low- 
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Fig.  5.  Trajectory  of  drifter  J.tJ4  (large  arrow)  and  SST  data  (  C)  for  June  29  to  July  5. 1915.  Shorter  arrows  denote 
the  flow  at  the  edges  of  the  Loop  Currerit  while  dotted  tines  indicate  the  locations  of  Loop  Current  waters  (tased  on  XBT 
data). 
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Fig.  6.  Ill)  Cruise  tracks  for  XBT  daia  collected  during  July  Il>-I9.  1985.  The  arrows  denote  the  flow  at  the  edges  of 
the  l.oop  Current  and  ring  based  on  the  vertical  temperature  structure  .shown  in  (h)  and  (r)  both  temperature  data  from 
the  M/V  Ndi  Co  6  cruise,  July  16-19.  1985.  and  (d)  temperature  data  from  the  5rcmi  Hispanig  cruise,  July  16-17, 1985. 


pass  liltcr  in  making  their  calculation t  for  a  ring  in  the  west¬ 
ern  GOM.  it  was  found  in  this  study  that  such  a  liltcr  was  not 
sufficient  to  remove  higher-frequency  fluctuations  recorded  by 
drifter  3.754.  As  a  result,  we  used  a  164-hour  (half  power  point) 
low-pass  filter  on  the  velocity  data  of  the  drifters.  We  then 
calculated  the  period  of  rotation,  the  ellipticity,  the  orienta¬ 
tion.  am’  the  velocity  about  the  translation  ring  center  (swirl 
velocity)  as  seen  by  drifters  3354  and  .3378. 

Ijiop  Current  Kinematics 

The  observed  (filtered)  Loop  Current  velocities  as  seen  by 
drifter  3354  are  shown  in  Figure  12.  The  speed  fluctuates  in 
magnitude  from  15  to  80  cm/s,  and  there  are  minimums  at  the 


beginning  of  tiie  record  and  during  Julian  days  210-225.  It 
was  during  days  220-230  that  drifter  3354  moved  northwest¬ 
ward.  aiT  apparent  rapid  northward  extension  of  the  Loop 
Current.  The  swirl  velocity  associated  with  the  rotation  of  the 
Loop  Current  is  shown  in  Figure  13.  These  magnitudes  also 
vary  from  1 5  to  80  cm/s. 

The  extent  of  water  involved  with  the  northward  movement 
of  the  Loop  Current  is  indicated  by  Figure  14,  which  shows 
the  time  history  of  the  distances  from  the  drifter  to  the  center 
of  rotation.  Prior  to  the  northwest  movement,  the  maximum 
radius  of  the  rotation  as  seen  by  the  drifter  was  8S  km.  After 
the  northwest  movement  the  radius  increased  up  to  100  km. 

The  rotational  frequency  of  the  Loop  Current  is  shown  in 
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Fig.  7.  Tfajeciofie^  of  drifters  i354  fsolid  targe  arroui  and  3378  (dashed  large  arrow)  and  SST  data  (  O  for  July  16-23. 
IIXS.  Sluiricr  arrows  dcnoic  the  (low  at  the  edges  of  the  l.oop  Current,  while  dotted  lines  indicate  the  locations  of  Lood 
Current  waters  (based  on  .XBT  data). 


Figure  15.  (For  the  calcul.ition  of  this  term,  .sec  Kirwttn  ci  til.  ino\ed  ttorthwesl.  the  mean  period  of  rotation  decreased  from 
ri9K4].)  Figure  15  .show.s  a  general  decrease  and  then  increase  about  1.3  day.s  to  approximately  lO.S  days.  If  we  compare 
in  frequency,  with  a  minimum  occurring  right  before  the  Figures  U  and  15.  we  see  that  longer  periods  tend  to  be 
^•lOtthwt^■twal<i  extension  of  the  water  mass.  As  the  water  associated  with  larger  radii. 
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Fig.  ».  1(1)  Cruise  track  for  XBT  data  ctillectcd  during  August  16-18.  1985.  The  arrow  denotes  the  (low  at  the  edge  of 
the  Loop  Current  based  on  the  vertical  temperature  structure  shown  in  (61  lemperature  data  from  the  Sieiu  llispania 
cruise.  August  16-18. 198.^ 
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I'ig.  9,  Trajectories  of  drifters  3.154  (solid  large  arrow)  and  3378  (dashed  large  arrow)  and  SST  data  ('C)  for  August 
12  19.  (985.  .Shorter  sirrow  denotes  the  How  at  the  edge  of  the  Loop  Current,  while. dotted  lines  indicate  the  location  of 
Loo|i  Current  waters  (based  on  ,XDT  data). 

As  for  the  .shape  of  the  lielti  of  rotation,  we  tise  eccentricity  Kincimtiics  ofRinu  3S7ft 
c.  defined  as  the  major  a\is  length  divided  by  the  minor  axis 

length.  Also,  we  calctilittc  the  north.'cast  orientation  of  the  The  filtered  velocities  of  ring  3.378  are  shown  in  Figure  17. 
major  axis,  and  both  of  thc.se  vtiriablcs  arc  .shown  in  IHgure  The  magnitudes  of  these  oscillations  are  relatively  large,  up  to 
16.  The  Moss  field  starts  out  rather  elliptical  (<•  f  1.7.5)  but  ~  85  cm/s  with  a  minimum  of  about  50  cm/s.  The  swirl  spieds 
becomes  more  circular  by  day  250  (c  =  1.4).  The  elliptical  are  shown  in  Figure  18.  and  we  see  here  an  initial  increase 
orientation  was  mostly  cast  west  at  first  but  eventually  from  20  to  75  crit/s  followed  by  a  decrease  to  about  50  to  60 
became  ntorc  northwest  southeast.  cm/s.  These,  variations  in  swirl  magnitude  coincide  with  vari- 


Hg  10  Trajectories  of  drifters  33.S4  (solid  large  arrow)  and  3378  (dashed  large  arrow)  and  SST  data  ('Q  for  September 

10-17.1985. 
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Fig.  1 1.  TrajccloriM  of  drifters  .1.154  (solid  large  arrow)  and  3378  (dashed  large  arrow)  and  SST  data  CQ  for  September 

17-24.  1985. 


iitioru  in  the  size  of  the  circle  of  rotation  (Figure  I).  Figure  19 
quahtilies  the  distances  from  the  center  of  rotation,  with  a 
ma.simum  radius  0195  km. 

The  rotation  frequency  of  ring  3378  is  shown  in  Figure  30. 
1  he  initial  minimum  in  frequency  ( -0.1 1  ' ')  is  followe<l,''v  an 
increase  to  about  0.14"'.  From  about  day  220  and  on,  the 
rotational  frequency  of  ring  3378  gradually  decrexsed.  with  (he 
period  of  rotation  increasing  from  about  8  days  to  about  9 
days. 


t.oop  Curren*  (3354) 


Julian  Days  (1985) 

Fig.  12.  Time  histories  of  the  nitered  observed  speed  components 
for  drifter  .1354.  The  light  curve  is  the  easi/wesi  speed,  and  the  darker 
curve  is  the  nsirih/soulh  .speed. 


Con.sidering.ihe  shape,  the  data  indicate  that  the  ring  was 
initially  rather  elliptical  (Figure  21).  Within  two  rotation 
periods  (his  eccentricity  decreased  somewhat,  with  (he  ring 
becoming  almost  circular.  The  results  shown  in  Figure  21 
imply  that  the  initial  ellipse  orientation  was  mostly  east/west. 
This  orientation  did  not  seem  to  change  significantly  as  the 
ring  became  more  circular. 

5.  Discussion 

The  cliissiciil  concept  of  Loop  Current  processes  is  one  in 
which  part  of  the  Guif  Stream  at  first  flows  directly  from  the 
Yucatan  Straits  to  the  Florida  Straits.  Within  this  flow  field. 


I'ip.  1.1.  Time  hi.viories  of  the  swirl  speed  components  about  the 
center  of  roiation  for  ihe  Loop  Current.  The  light  curve  is  the  east/ 
west  sivccd.  :ind  (he  d.irkcr  curve  is  Ihe  norlh/souih  speed. 
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Fig.  14.  Time  histories  of  Ihe  east/west  (light  curve)  and  north/ 
south  (darker  curse)  distances  of  drifter  3.'54  from  the  center  of  rota¬ 
tion  in  the  Loop  Current. 
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Fig.  16.  Time  histories  of  the  Loop  Current  eccentricity 
(asterisks)  and  orientation  of  major  axis  (dashed  lines).  For  the  flow 
fleld  orientation,  north  is  in  the  positive  t  direction,  and  east  is  in  the 
positive  time  direction. 


inMabilitic.s  exist  which  result  in  meandering  of  the  Loop  Cur¬ 
rent  llliirlbiirt  ami  Thompson,  1980].  As  the  sues  of  the  mean¬ 
ders  increase  and  reach  northward,  it  is  generally  believed  that 
the  now  Held  wraps  back  onto  itself  and  “shorts"  across  the 
stream  of  (low:  part  of  the  How  would  still  go  northward 
around  the  Loop  Current  extension,  while  the  remainder  of 
the  now  would  take  the  more  direct,  southerly  route  to  the 
Florida  Straits.  This  is  analogous  to  other  geophysical  phei 
nomena.  such  as  the  creation  of  ox-bow  lakes  by  meandering 
rivers.  I■■inally.  as  the  curvature  of  the  How  reaches  its  maxi- 
tmiiii.  iiioic  of  the  l.oup  Cuirciil  Mow  lakes  the  soulherly 
route,  and  a  OOM  ring  is  eventually  pinched  olT. 

The  interesting  point  indicated  by  the  data  prc.sehled  here  is 
the  closed  anticyclonic  How  within  the  Loop  Current  immedi¬ 


Fig  15.  Time  history  of  the  rolalioiial  frequency  within  the  Loop 
Current. 


ately  after  ring  3.^78  was  pinched  olT.  Also,  the  trajectory  of 
drifter  3354  reveals  (hat  this  rotational  feature  initially  existed 
oT  the  coast  of  Cuba.  (I  has  been  shown  in  numerical  studies 
that  the  How  of  the  Loop  Current  can  itself  create  an  anti¬ 
cyclonic  How  Heid  olT  (he  northwestern  coast  of  Cuba 
[Thompson,  1986].  The  northward  flow  west  of  the  (ip  of 
Cuba  along  with  the  southeasterly  flow  along  the  north  cen¬ 
tral  coast  of  Cuba  will  obviously  produce  negative  vorticity  ofl* 
Cuba’s  northwestern  shore.  Such  an  anticyclonic  flow  has 
Ix'cn  documented  by  a  number  of  hydrographic  surveys  [e.g., 
Nowlin  mul  McLellan,  1967;  Cochrane,  1972].  Using  hydro- 
graphic  data.  Hofmann  and  K'or/ey  [1986]  estimated  Ihe  depth 


Fig.  17.  Time  histories  of  the  filtered  observed  speed  components 
for  ring  3.378.  The  light  curve  is  the  easi/west  speed,  and  the  darker 
curve  is  the  north/south  speed. 
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I  iiL  IS.  Time  hi<h>rics  ••!  «lw  «i»ccU  t5*fiipJ=snis  aKml  the 
center  t<(  rKlnlion  for  riiic  .v*7S  'I  he  llchi  cane  i*  the  esM  »crt  '>fccil. 
:iii(J  the  (l.-irLer  cnric  i«  the  north  <<>uth  epcni. 

t«r  ihi':  Cithttn  eddy  to  he  -^kIXI  hi  uiih  a  current  maenitude  ttf 
IO-25  cm  c.  in  good  agreement  with  the  data  frcim  this  study 
(l-iutirc  12).  Also,  data  .show  that  lltc  temperature  and saiinili* 
characteristics  of  GOM  rings  anti  the  Cuban  eddy  are  practi¬ 
cally  identicjil  rAfi7.t//«i»i.  IWirt:  Hllini.  19821. 

I  he  Cuban  eddy  appears  to  be  the  flow  lield  in  which  drif¬ 
ter  .l.l?4  initially  became  entrained.  'Flic  drifter  trajectory  indi- 
ctties  that  this  Ilow  lield  moved  slightly  northwest  while  ring 
.t.t  iX  was  still  connected  to  the  I.otip  Current.  Hut  it  was  only 
slioriiy  after  ring  .7.7 /X  pinched  olT  that  this  anticyclone  pushed 
more  :ii»rthwe<lerly.  to  about  26,.s  N.  Figures  10  and  1 1  indi¬ 
cate  that  after  this  ftnal  push  (late  August  19X5).  the  rotational 
liclil  bccauic  ail  integral  component  of  the  Loop  Current, 
riiiis  the  data  in  (his  c.ise  imply  that  the  kernels  of  Loop 
('iirrent  rings  may  come  from  waters  olT  the  northwest  coast 
of  Cuba. 


Julian  Days  mfl'i) 


I  ig.  I')  Time  histories  of  the  east  west  ( ighi  ciirscl  and  north/ 
votiih  Idarkcr  curve)  distances  of  drifter  .t.778  Irorii  the  center  of  rota¬ 
tion  of  ring  .177X, 


Julian  Cavs  119*51 

fig.  2<Jl  'Tmiclttitory  of  die  rotational  frequency  of  ring  2378. 

An  e.sample  of  a  modeled  ring  pinch-oAT  in  the  GOM  is 
shown  in  Figures  22  and  23  (from  Wtrikri^i,  1986].  On  day 
63  (Figure  22)  tlie  model  shows  a  configuration  similar  to  that 
when  drifter  .7354  ants  seeded.  The  ting  is  well  defined  as  the 
northern  e.sfcnsion  of  the  Loop  Current,  but  it  has  not  com¬ 
pletely  pinched  olT.  Note  that  there  docs  exist  a  hint  of  anti- 
cyclonic  motion  off  the  northwestern  coast  of  Cuba.  By  day 
153  f3  months  later)  the  model  shows  the  ring  free  of  the  Loop 
Current.  Morcovtr.  the  anticyclone  off  Cuba  is  larger  and 
much  better  defined,  ilowiever.  the  modd  shows  the  northern 
.edge  of  the  Loop  Current  extending  to  only  26'N.  whereas 
XBT  and  drifter  data  presented  here  show  an  extension  to 
27  N  in  a  .simiLir  lime  period  (3  months). 

It  IS  not  until  after  9  months  before  the  Loop  Current  in 
this  .simulation  pushes  northward  again  in  the  process  of 


Julian  Days  (1985) 

iTg.  21.  'lime  li!siorie.s  of  the  eccentricity  (asierbk.s)  and  major 
.'ISIS  oricni.'iliiin  (dashed  line.s|  for  ring  3378.  For  the  flow  field  orien- 
laiiun.  north  is  in  the  po.siiive  e  direction,  and  ea.si  is  in  the  [>ositive 
lime  dirccliiin. 
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liiiicliiiig  oir  aimllwr  ling  (l-igiire  23».  By  ihis  lime,  ihc  aiili- 
ctckinc  ihal  hail  heen  spun  un  by  Ihc  Lmip  Currcnl  ha^  Io*l 
its  kicniily  ic!:iy  7??.  I'iijiirc  2.7|.  Thiw  Ihc  priKCS*  a<  dcler- 
r.:ii:c:l  fi-Tn  !!  ii-:=<'ry  I'f  <::;f'cr  .■*.754  is  ^l•l^■crcn^  from  Ihis 
;-:;r:icu-:;>‘  <'ii  iwo  major  ways.  I'irsl.  Ihc  Loop  Cor- 

...  vcj'  •■;i.ns:  i!iiiiliw-;,.'ti  lo  27  N  in  2.5  months 

'..I  ‘::!av  .M  10  r-H'i!thS  hlf  Mldl.niorfcl 
M.-::..  .Seci*:i-,-.  !!•!■:  rapM  -.•.siCKsion  was  <iccbm|miiicd  by 
circii!a"0!!  voiict  iirigiiiallv  gcncralcd  oir  Ihc  coast 
.1  5. 'ah...  JliAvovci.  li!  other  simulations  with  no  d«p  in-flow 
ihiough  ilic  Yucatan  .Straits,  signilicaiU  Loop  Current  pene¬ 
trations  can  iHX'iir  in  40  days.  Clearly,  this  is  a  phenontiena 
that  requires  Turlhcr  study. 

CurrvK!  Vivxtix  Rhitf  Khicnuahs 

I^rifters  .1354  and  3378  were  both  in  the  GOM  during 
Julian  days  199  264.  1985.  However,  the  initial  variations  of 
the  kinematic  parameters  of  ring  3378  (Figures  17-211  all  indi¬ 
cate  an  adjustment  periinl  from  days  l‘W  220  during  which 
lime  the  drifter  (with  its  2(IO-m  wciglilcrl  linci  found  an  appro¬ 
priate  orbit  williin  the  ring.  Thus  we  compare  the  ring  and 
Loop  Current  kinematics  from  day  225  trr  day  264.  At  the 
ix’ginning  of  tin':,  time  period  Ihc  Loop  Current  was  rotating 
wiih  I  period  of  aSoul  10.5  days,  while  ring  3378  had  a  perioil 
of  S.2  liii.v:;.  The  oricniations  of  the  ellipses  of  rotation  ll•■ig- 
;:ra.s  lit  a'ul  21)  were  quilesimilar  igcnerally  cast,  westl.  butThe 
L.v'j'  Curriiii  was  slightly  more  elliptical  (c  »-  L5|  than  ring 
3.'7?'- !<■*  1.251. 

I  he  .vwirl  speeifs  am!  the  'radii  of  revolution  for  the  drifters 
in  hoih  aniicychmcs  were  quite  similar  for  days  225-2(4:  80 
cm  s  l.'r  the  swirl  spccrls  ami  '  9(!  km  for  the  r.idii  relative  to 
the  rii'.,!  center,  luiereslingly.  both  aiilii;!. clones  e.ipeiienced 
lo'nao."  yetiods  of  rolaiidn  lor  smaller  disj.tiices  from  the  center 

,.f  iitiiiJioa. 

Kiiu-niiiHiX  oiul  i/rcc  a  Riii<i  Oc/uc/.nicnr 

Wv  now  iiirii  lo  the  kinematics  of  the  Loop  Current  at  its 
mos!  iKuiherly  e.xicnsion  (Julian  days  2,30-264)  and  the  kin- 
-.•laa-iic'  of  ring  3378  .soon  after  it  broke  free  from  the  Lot>p 
Cm  rent  (Julian  days  225-235).  Keeping  m  mind  that  thew 
data  represent  two  difl'erent  anticyclonic  phenomena,  how 
typical  arc  the  COM  rings  before  and  after  they  break  off 
from  Ihc  Loop  Current? 

The  similarities  between  Ihc  two  sets  of  kinematics  arc  con¬ 
siderable.  First,  the  magnitudes  of  the  swirl  speeds  and  of  the 
radii  arc  practically  identical.  Moreover,  their  elliptical  orien¬ 
tations  arc  both  norlhwcsi/southcasl.  the  Loop  Current  was 
slightly  more  ciliplibt  than  ring  3378.  but  one  might  expect 
such  a  dill'crcncc  in  ciliplicily  seeing  that  an  anticyclone  still 
.'iiachcd  lo  Ihc  Loop  Current  is  a  "capluicd”  phenomena, 
while  a  detached  ting  is  an  isolated  vortex.  The  similarities  of 
the  anticyclones  plus  their  close  periods  of  rotation  imply  that 
the  basic  kinematic  characteristics  of  GOM  rings  can  be  es¬ 
tablished  as  the  Loop  Current  pushes  northwestward  off  Ihc 
shore  of  Cuba. 

Smimuiry 

riie  scenario  implied  here  is  one  in  which  the  essence  of  the 
ring  is  established  ofl'the  northwest  coast  of  Cuba,  develops  in 
deeper  waters  olV  the  coastline,  and  then  lakes  on  its  final 
characteristics  shortly  after  the  previous  ring  coiuplelely  de¬ 
taches  from  the  Loop  Current. 


The  implicaiion-of  .such  a.  mechanism  h  .that  of  time  scales. 
One  is  not  required  to  wail  until  instabilities  in  llie  flow  field 
grow  large  enough  to.produce  a  closed  rotational  feature.  The 
data  heie  show  that  a  new'  Loftp  ('uireiit.rotalionai  feature 
can  lx  well  esiabiishesi  even  before  the  previous  ring  has  to¬ 
tally  pinched  ofl'.  It  was  established  that  ring  3378  was  de¬ 
tached  Ironi  the  Loop  Curfen:  by  iuid-July  l!/85.  Yet  it  wiis 
seen  that  the  new  rotational  structure  had  moved  to  latitude 
26.S'N  by  the  beginning  of  Scpiemher  1985.  thus  the  Loop 
Current  wa.s  reconfigured  for  another  ring  separatidn  only  1.5 
months  after  the  previous  ring  separation.  Elliqi  [1982]  docu- 
mcAled  three  ring  separations  in  a  12-monlh  period  using 
hydrographic  data.  The  presence  of  Ghost  Eddy  along  with 
ring  3378  implies  a  similar  sc|Kiralion  rate.  Model  studies 
(c.g.,  Hufihuri  and  Thompson,  1980;  Thompson,  1986;  Wall- 
m(fl>  1986]  report  separation  rales  ranging  from  4  to  18 
months,  depending  upon  the  type  of  instability  involved. 
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Observed  and  Simulated  Kinematic  Properties  of  Loop  Cunrent  ^ngs 

A.  b.  Kirwan,  J.  K.  Lewis,^  A.  W.  Indest*  -P.  Reinersman/  and  I;  Quintero* 

Two  rings,  shed  by  the  Loop  Cuneot  in  1980  and  19824-  were. observed' for  several  montto'by 
satellite-tracked  drifters  to  migrate  acro»  the  Gulf  of  Mexico.  The  drifter  path  data-have  been  invert^ 
to  obtain  estimates  of  the  paths  of  the  centers  of  the  two  rings,  ring  shape,  and  the  swri  velocities.  Three 
drifters  were  deployed  in  the  1980  ring,  and  the  analysis  of  that  data  set  establishes  the* variability  of  the 
above  kinematic  estimates  for  one  ring.  A  comparison  of  the  analysis  of  data  from  both  rings  provides 
some  idea  on  inter-ring  variability.  Both  rings  impacted  the  Mexican  continental  slope  at  about  218°N, 

95.S’W.  After  a  brirf  adjustment  period,  both  rings  reestablished  and  maintained  a  vortex  character  for 
several  months  in  the  slope  region  while  migrating  slowly  to  the  north.  The  paths  of  the  centers  of  the 
two  rings  along  the  slope  are  virtually  identical.  The  stme  analysis  routine  was  applied  to  some  simulat¬ 
ed  drifter  data  obtained  from  the  Hurlburt  and  Thompson  (1980)  Gulf  of  Mexico  primitive  equation 
model.  In  the  midgulf,  the  agreement  between  the  observed  rings  and  the  simulate  ring  is  good, 
although  the  fonher  showed  stronger  interaction  with  the  continental  slope  topography  and/or  circu¬ 
lation  than  was  seen  in  the  latter.  Along  the  slope,  the  model  ring  kinematic  char^eristics  were  in 
extraordinary  agreement  with  the  observations. 


1.  Introduction 

The  circulation  of  the  Gulf  of  Mexico  has  been  the  subject 
of  a  surprising  amount  of  theoretical  and  observational  stud¬ 
ies.  From  the  turn  of  the  century  [Sweitzer,  1898]  until  1973 
[Austin,  l95S’,.Nowlin  and  McClellan,  1967;  Nowlin,  1972],  the 
emphasis  was  on  analyses  of  hydrographic  data.  These  studies 
established  the  presence  of  poorly  defined,  but  quui- 
permanent,  anticyclonic  signatures  in  the  central  part  of  the 
gulf. 

Since  1973,  the  emphasis  has  shifted  from  qualitative  de¬ 
scriptions  of  the  hydrography  to  attempts  at  quantifying  dy¬ 
namical  mechanisms  associated  with,  the  anticyclonic  struc¬ 
tures.  Sturges  and  Blaha  [1976]  and  Blaha  and  Sturges  [1978] 
proposed  that  the  western  gulf  circulation  was  the  result  of  the 
balance  of  wind  stress  curl  and  planetary  vorticity,  i.e.,  a  mini¬ 
western  boundary  current  effect  in  the  Gulf  of  Mexico.  Elliot 
[1979,  1982]  disputed  this,  noting  that  this  balance  was  not 
generally  achieved  if  synoptic  wind  data  on  a  fine  scale  were 
used.  His  work  also  gave  more  credence  to  the  earlier  idea  of 
Ichiye  [1962]  that  the  anticyclonic  features  were  in  fact  rings 
that  had  been  shed  by  the  Loop  Current  in  the  eastern  gulf. 
Kirwan  et  al.  [1984u]  substantiated  this  hypothesis  by  track¬ 
ing  the  migration  of  a  ring  shed  by  the  Loop  Current  across 
the  gulf  to  the  continental  shelf  off  of  Mexico. 

In  a  companion  paper,  Kirwan  et  al.  [19846]  assessed  the 
translation  velocities,  local  vorticity,  deformation  rates,  and 
shape  of  the  ring  as  it  propagated  across  the  gulf.  The  analyses 
of  the  three  drifters  that  were  in  the  ring  simultaneously 
showed  generally  good  agreement  of  these  properties.  The  re¬ 
sults  also  are  in  general  qualitative  agreement  with  the  nu- 
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merical  simulations  from  the  general  circulation  model  (GCM) 
of  Hurlburt  and  Thompson  [1980],  hereinafter  referred  to  as 
HT. 

Is  the  qualitative  agreement  between  GCM  calculations  and/ 
one  particular  ring  fortuitous  or  the  result  of  deterministic 
dynamics?  We  examine  this  question  with  observations  from 
two  different  rings  along  with  simulated  data  from  the  HT 
model.  Specifically,  we  compare  in  detail  kinematic  properties 
such  as  the  movement  of  ring  centers,  ring  translation  and 
swirl  velocities,  and  ring  shapes  from  198(3-1981,  1982-1983, 
and  simulations  from  the  HT  GCM. 

The  analysis  routine  used  to  determine  these  kinematic 
properties  is  an  improved  version,  of  that  reported  by  Kirwan 
et  al.  [19846].  The  routine  inverts  L)i;;<angian  path  data  to 
obtain  the  desired  kinematic  properties.  Unlike  the  earlier 
study,  the  present  routine  emphasizes  ring  shape  and  trans¬ 
lation.  Also,  considerably  better  time  resolution  on  parameter 
estimates  can  be  made  with  the  new  routine.  Some  details  are 
given  below  and  in  the  appendix.  The  internal  consistency  of 
the  analysis  routine  is  being  assessed  separately  with  simulat¬ 
ed  data  from  the  HT  GCM. 

The  following  observed  data  were  utilized  in  this  study. 
These  come  from  two  rings  which  occurred  in  1980-1981  and 
1982-1983.  We  have  reported  previously  on  the  former  [see 
Kirwan  et  al.,  1984o,  6]  using  the  more  primitive  analysis 
routine.  That  ring  had  three  drifters  in  it  simultaneously,  and 
so  our  re-analysis  provides  some  indication  of  intraring  varia¬ 
bility  along  with  some  indication  of  the  consistency  of  the 
analysis  routine.  Only  a  descriptive  analysis  of  the  latter  ring, 
which  had  only  one  drifter,  has  been  reported  before  [Lewis 
and  Kirwan,  1985].  Comparison  of  these  two  rings  provides 
some  measure  of  inter-ring  variability. 

These  observed  data  are  compared  with  simulated  La- 
grangian  data  from  the  HT  GCM.  This  model  was  a  two-layer 
(lower  layer  active)  nonlinear  primitive  equation  model  of  the 
Gulf  of  Mexico  with  realistic  bottom  topography.  The  grid 
spacing  for  the  calculations  used  here  was  10  km.  The  simulat¬ 
ed  path  was  obtained  by  advecting  a  parcel  at  each  model 
time  step  (4  hours)  to  a  new  position  given  by  the  product  of 
the  velocity  vector  at  the  current  position  of  the  parcel  and  the 
time  step.  Since  the  position  of  the  parcel  was  rarely  at  a  grid 
point,  it  was  usually  necessary  to  interpolate  the  velocity 
vector  from  surrounding  grid  points.  These  calculations  were 
kindly  performed  by  A.  Wallcr^t. 
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Both  the  observations  and  the  simulations  track  appropri¬ 
ate  rings  for  many  months.  During  this  time,  the  real  and 
simulated  rings  move  from  the  midgulf  region  with  fairly  flat 
topography  in  excess  of  3000  m  to  the  continental  slope  off  of 
Mexico.  This  area  of  interest  is  depicted  in  Figure  1.  In  the 
midgulf  region,  ring  translations  are  governed  essentially  by 
Rossby  wave  dynamics,  and  so  one  expects  them  to  move 
west-southwest.  But  along  the  slope,  interactions  of  the  ring 
with  bottom  topography  become  important.  Since  the  dynam¬ 
ics  in  these  two  geographic  regions  are  so  different,  we  have 
performed  separate  analyses  for  the  rings  in  the  two  regions. 


Fig.  1.  Depiction  of  the  area  in  which  the  observed  and  simulated 
Loop  Current  rings  are  studied. 


2.  Analysis  Proced'Ire 

The  observations  used  in  this  study  are  the  positions  of  the 
Polar  Research  Laboratory  drifters  with  Argos  identifications 
1598,  1599,  1600,  and  3374.  Kirwan  et  al.  [1984<j,  b]  reported 
on  the  analysis  of  1598,  1599,  and  1600,  which  were  deployed 
simultaneously  in  November  1980.  Position  data  from  these 
three  drifters  are  analyzed  to  establish  the  variability  within  a 
ring.  Intraring  variability  is  assessed  by  comparing  the  analy¬ 
sis  performed  on  the  three  drifters  in  the  1980-1981  ring  with 
that  of  drifter  3374,  which  was  deployed  in  the  1982-1983 
ring.  Each  drifter  was  the  standard  Polar  Research  Labora¬ 
tory  buoy  drogued  by  a  200-m  thermistor  line.  No  thermistor 
data  are  available  from  any  of  these  units. 

The  analysis  procedure  utilizes  a  parametric  kinematic 
model  discussed  by  Kirwan  et  al.  [1984b]  for  inferring  the 
translation  and  swirl  velocities  as  well  as  the  vorticity,  hori¬ 
zontal  divergence,  and  deformation  rates.  This  model  was  first 
proposed  by  Okubo  [1970]  to  describe  small-scale  motion 
about  flow  singularities.  Kirwan  [1984]  showed  that  same 
model  can  be  applied  to  larger  scales  if  the  model  parameters 
arc  constant  along  a  parcel  path. 

The  model  assumes  that  the  swirl  velocity  is  given  as  a 
linear  function  of  the  distance  from  a  local  flow  center.  Here 
the  local  parcel  flow  center  will  be  referred  to  as  the  “ring 
center”  or  soipe  times  “orbit-determined  center”.  This  termin¬ 
ology  is  dictated  only  by  the  nature  of  the  available  data. 


Flierl  [1981]  has  pointed  out  that  this  center  may  not  coincide 
with  an  “Eulerian”  center.  (We  interpret  this  as  a  center  of 
mass.)  Smith  and  Reid  [1982]  have  developed  a  general  for¬ 
malism  for  studying  the  centers  of  the  distribution  on  any 
property,  for  example,  mass,  potential  energy,  kinetic  energy 
or  enstrophy.  In  general,  these  centers  do  not  coincide,  nor  do 
their  propagation  velocities,  ^tablishing  the  relations  be¬ 
tween  these  various  centers  for  Gulf  of  Mexico  rings  requires 
much  more  data  than  is  now  available. 

For  the  present  model  the  total  velocity  (u,  u)  of  a  drifter  is 
composed  of  the  translation  velocity  (Uj-.  Vj)  of  the  ring 
center  and  the  swirl  velocity  (Uj,  Uj): 

U=Vy  +  Us 
v=Vt  + Vs 

Uj  =  (d  +  a)x/2  +  (b  -  c)y/2 

(2) 

Uj  =  (b  c)x/2  +  (d  -  a)y/2 

Here  x  and  y  are  the  instantaneous  coordinates,  relative  to  the 
ring  center  of  a  particular  parcel  (drifter).  If  (2)  is  regarded  as  a 
Taylor  expansion  for  the  velocity  field,  then  the  parameters  a 
and  b  can  be  interpreted  as  deformation  rates  and  c  and  d  as 
the  vertical  vorticity  and  horizontal  divergence.  It  is  stressed 
that  these  parameters  are  evaluated  along  individual  paths. 
Eulerian  estimates  of  deformation  rates,  vorticity,  and  diver¬ 
gence  likely  would  differ  from  these  estimates.  Consequently,  it 
is  better  to  interpret  a,  b,  and  c  as  shape  functions  which 
determine  the  orientation  and  ellipticity  of  a  particular  orbit 
in  a  ring.  See  Kirwan  et  al.  [1984b]  for  details. 

The  general  solution  to  (1)  and  (2)  is  easily  obtained 
lOkubo,  1970].  The  procedure  here  is  to  apply  this  solution  to 
every  time  interval  between  fixes  (in  practice,  interpolated 
positions).  Thus  for  the  interval  t,i^  t  £  t^^.^,  one  obtains 

“*  =  +  (<^k  +  «k )V2  +  (6k  -  (3) 

fk  =  +  (6k  +  c^)x^/2  +  (dj  -  aM2  (4) 


where 

Xk  =  {exp  [ru(»-«k)3[2fk(y‘»k^  +  6k^-Ct*+a»)+  yk(6k-Ci)] 
-exp  [r«(t-tj)][2ft(at- Va7+V^^)+  yk(6k-Ck)]} 

2Va*^+6**-c»^  (5) 

^k  =  {exp[r,»(r-g][;fA+Ck)+  nOk^  +  6k^-Ck^-Uk)3 

+exp  [r2*(f-t»)][-2f»(6*+c*)+  Tk(\/V+V^^+«k)]} 

-  2v/ak^  +  6**  -  c»^  (6) 


In  (5)  and  (6),  (A’t,  are  the  coordinates  of  the  parcel  in 
question  at  time  t  -  t*  relative  to  the  ring  center.  Also, 

r,k  =  (d*  +  +  6*^  -  c^^)/2  (7) 

'•a  =  (‘'k-v/V  +  6»^-c*^)/2  (8) 

are  the  eigenvalues  of  the  matrix 


2M»  = 


■(dk  +  Ok) 
.(6k  +  Ck) 


(6k-Ck)‘ 
(dk  ^  Ok). 


Note  that  these  eigenvalues  are  complex  conjugates  when 
Ck*  >  +  b^^.  This,  naturally,  produces  real  periodic  solu- 
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tions  which  are  a  characteristic  of  swirl  velocities  of  drifters 
trapped  in  a  ring. 

To  apply  this  model  to  drifter  data,  differentiate  (3)  and  (4) 
three  times  with  respect  to  time.  One  side  of  these  equations 
contains  the  unknown  variable  set  =  (a*,  h*.  c*.  d*.  .Y*.  X). 
The  other  side  of  these  equations  can  be  estimated  from  path 
data.  These  six  equations  do  not  contain  the  translation  veloc¬ 
ity,  since  it  is  assumed  a  priori  to  be  constant  for  a  time 
interval.  (But,  of  course,  its  value  may  change  from  one  time 
interval  to  the  next.)  These  six  equations  are  inverted  analyti¬ 
cally  for  the  six  elements  of  Z^.  The  translation  velocities  can 
be  recovered  from  (3)  and  (4),  since  all  other  quantities  are 
now  known.  Alternatively,  the  translation  velocities  can  be 
obtained  by  differentiating  the  path  of  the  ring  center.  Both 
procedures  provide  comparable  estimates.  Here  we  use  the 
latter  procedure.  With  all  elements  of  Z^  now  determined,  one 
then  moves  to  the  next  time  interval  and  repeats  the  calcula¬ 
tions,  thus  obtaining  a  time  series  of  Z^.  Details  are  given  in 
the  appendix. 

We  have  expended  a  great  deal  of  effort  in  learning  how  to 
make  reliable  time  derivative  estimates  of  the  path  data.  For 
the  data  analyzed  here,  the  raw  position  files  have  been  edited, 
interpolated  to  equally  spaced  time  intervals,  and  low-pass 
filtered  (half  power  point  of  100  hours).  All  derivative  esti¬ 
mates  were  obtained  by  fourth-order  accurate,  centered  differ¬ 
ence  schemes.  Each  derivative  file  was  edited  and  low-pass 
filtered.  The  calculations,  as  described  in  the  appendix  for  Z^, 
were  performed  with  the  latter  derivative  files.  The  calculated 
Z^,  in  turn,  were  edited  and  low-pass  filtered. 

The  analysis  procedure  just  outlined  is  an  improvement 
over  that  used  by  Kirwan  et  al.  [19846].  Their  method  re¬ 
quired  an  a  priori  assignment  of  a  solution  form  to  (4)  and  (5); 
i.e.,  (7)  and  (8)  were  assumed  to  be  complex  conjugates,  and 
values  of  Z^  were  assumed  to  be  constant  over  one  ring  revo¬ 
lution  (approximately  15  days).  The  procedure  used  here  does 
not  require  an  a  priori  assignment  of  the  solution  form;  it  lets 
the  data  determine  this.  Thus  it  should  have  wider  utility.  In 
addition,  the  Z^  are  now  required  to  be  constant  over  seven 
time  intervals  (42  hours  in  the  present  case)  rather  than  2 
weeks. 

As  was  indicated  above,  a,  b,  and  c  are  not  especially  useful 
for  comparison  with  Eulerian  data.  For  rings  a  more  appro¬ 
priate  description  would  be  the  elliptical  structure  of  the  orbit 
traversed  by  a  particular  drifter.  As  was  shown  by  Kirwan  et 
al.  [19846j,  the  characteristic  equation  is 

(Cj  +  6j)xt^  +  (Cj  -  -  2ajXty» 

»  Lj  exp  [-d»(t  -  <»)]  (9) 

Here  is  the  angular  momentum  per  unit  mass  relative  to 
the  ring  center.  This  can  be  calculated  directly  (see  equation 
(A  18)).  In  (9)  the  argument  in  the  exponential  term  can  be 
made  zero  by  making  the  evaluation  at  the  beginning  of  the 
time  interval  t  =  t^.  Note  that  a  requirement  for  (9)  to  describe 
an  ellipse  is  that  -f  6*^. 

From  analytic  geometry  it  is  known  that  for  anticyclonic 
motion,  <  0,  the  major  axis  of  the  ellipse  makes  an  angle 
with  east  of 

aj  =■  i  tan"*  (-Ot/6j)  (10) 

In  evaluating  (10),  care  must  be  exercised  in  the  quadrant 
assignment.  Also,  the  semilength  of  the  major  axis  is 


(11) 

where 

The  paradigm  just  outlined  for  analysis  of  Lagrangian  data 
puts  critical  emphasis  on  time  derivatives  of  path  data.  Incon¬ 
sequential  errors  in  the  path  data  can  become  consequential  in 
the  derivative  estimates.  Moreover,  the  nonlinear  algebraic 
combinations  of  these  derivatives  in  the  paradigm  may  mag¬ 
nify  further  the  impact  of  these  errors  on  the  inversion.  The 
concern  is  that  two  paths  which  contain  the  same  kinematic 
information  but  differ  by  small  random  and/or  round  off 
errors,  will  produce  vastly  different  kinematic  estimates  when 
run  through  the  paradigm.  The  internal  consistency  of  esti¬ 
mates  of  ring  kinematics  is  presently  being  examined  with 
simulated  data  from  the  HT  GCM. 

3.  Midgulf 

3.1.  Observations 

Here  the  results  of  the  intraring  and  inter-ring  comparisons 
are  presented  for  the  midgulf  region.  This  includes  ring  kin¬ 
ematics  as  inferred  from  the  three  drifters  in  the  1980-1981 
ring  and  the  single  drifter  in  the  1982-1983  ring.  While  in  the 
midgulf  region,  both  rings  began  to  interact  with  the  conti¬ 
nental  slope  topography  and/or  circulation.  This  period  is  iso¬ 
lated  in  the  analysis.  The  same  analysis  routine  is  then  applied 
to  simulated  midgulf  ring  data  from  the  HT  GCM. 

The  1980-1981  ring.  The  data  for  this  portion  of  the  study 
come  from  three  drifters  (1598,  1599,  and  1600)  which  were 
seeded  in  the  1980-1981  ring.  Kirwan  et  al.  [1984fl,  6]  have 
already  reported  on  this.  The  data  have  been  reanalyzed  using 
the  general  algorithm  discussed  in  the  preceding  section  and 
in  the  appendix. 

Figure  2  shows  the  paths  of  these  three  drifters  in  the 
midgulf  region  along  with  the  positions  of  the  centers  and  the 
orbit  ellipses.  Arrows  on  the  drifter  and  center  paths  mark 
10-day  intervals,  while  the  ellipses  are  presented  at  15-day 
intervals.  The  paths  for  each  drifter  are  for  a  common  time 
interval  (day  340,  1980,  to  day  118,  1981).  This  figure  es¬ 
tablishes  three  points.  First,  all  three  make  the  same  swath 
through  this  portion  of  the  gulf,  suggesting  that  they  generally 
followed  the  ring.  The  overall  path  characteristics  of  drifters 
1598  and  1600  (Figures  2a  and  2c)  are  remarkably  similar, 
suggesting  that  they  were  in  very  similar  orbits.  Second,  the 
paths  for  the  center  of  the  ring,  as  determined  independently 
from  the  three  drifters,  is  in  good  agreement  as  far  as  about 
94°W  (the  first  2  months  of  the  record).  Discounting  the  very 
beginning  where  there  are  numerical  problems,  the  typical  dif¬ 
ference  in  the  location  of  the  center  at  any  one  time  is  30  km, 
about  20%  of  the  diameter  of  the  ring  as  determined  from 
satellite  imagery  [Kirwan  et  al.,  19846].  This  is  about  10  km 
larger  than  the  rms  displacement  between  the  maximum  sur¬ 
face  and  interface  displacement  anomalies  of  the  HT  GCM. 
The  30-km  variability  is  considerably  larger  than  that  report¬ 
ed  by  Hooker  and  Olson  [1984]  using  a  variant  of  the  tech¬ 
nique  employed  here  under  ideal  conditions.  In  view  of  our 
results  along  the  continental  slope,  described  below,  it  is  not 
clear  whether  the  large  variability  is  due  to  data,  technique,  or 
rapid  evolution  of  the  ring. 

Finally,  the  ellipses  also  show  generally  good  agreement 
east  of  94°W.  Note  the  tendency  to  develop  a  northwest  elon- 
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gation  as  the  ring  approaches  94°W.  West  of  94°W;  all  three 
orbits  suggest  that  the  ring  begins  to  intera«  with  the  conti¬ 
nental  slope  circulation  and/or  topography.  This  is  discussed 
later. 


LONGITUDE 
Fig.  2a 


LONGITUDE 
Fig.  2b 


Figure  3  shows  the  time  series  for  the  swirl  and  translation 
velocities.  All  three  records  indicate  swirl  velocities  in  excess  of 
50  cm/s  throughout  most. of  the  r^ord.  The  dominant  period 
is  about  14  days  with  a  variation  of  about  2  days  between  the 
records.  This  is  reco^ized  as  the  ring’s  rotation  period.  How¬ 
ever.  from  abburday  30  to  60.  1981.  all  three  drifters  show  a 
smaller-period  component  (about  5  days)  in  the  time  series 
and  a  decrease  in  swirl  amplitude.  As  will  be  discussed  below, 
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Fig.  2.  Drifter  paths,  inferred  center  path  (red)  and  orbit  ellipse 
(blue)  for  the  three  drifters  in  the  1980-1981  ring:  (a)  1598,  (6)  1599, 
and  (c)  1600.  Arrows  on  the  paths  are  at  104ay  intervals,  while  the 
ellipses  are  shown  at  15-day  intervals.  Here  S  and  F  refer  to  start  and 
finish. 
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Fig.  3.  Swirl  (left  scale)  and  translation  velocities  (nght  scale)  in¬ 
ferred  from  the  three  dnfters  in  the  1980-1981  nng.  (a)  1598,  (6)  1599, 
and  (c)  1600.  The  east  and  north  components  of  the  translation  veloc¬ 
ity  are  depicted  in  red  and  blue,  respectively. 
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it  is  speculated  that  during  this  period  the  drifters  were  re¬ 
sponding  to  effects  produced  by  the  continental  slope  and/or 
circulation  further  to  the  west.  All  the  translation  velocities 
show  a  consistent  westward  component  of  about  4  cm/s  up 
until  about  day  30,  1981.  For  the  reason  just  given,  it  is  not 
clear  that  the  analysis  applies  to  this  ring  for  the  period  from 
day  30.  1981.  to  day  60.  1981.  After  day  60,  1981,  the  ring 
centers  appear  to  be  stationary. 

The  1982-1983  ring  and  comparison  with  the  1980-1981 
ring.  The  same  characteristics  are  seen  in  the  analysis  for  the 
1982-1983  ring.  Figure  4  shows  the  drifter  path,  inferred 
center  path,  and  the  orbit  ellipses  for  drifter  3374.  As  in  Figure 
2,  there  is  significant  distortion  of  the  orbits  west  of  94°W. 
Unlike  the  earlier  ring,  no  significant  northwest  ellipticity  is 
developed  as  the  ring  approaches  94°W.  Comparison  with 
Figure  2  indicates  that  the  1982-1983  ring  center  followed 
nearly  the  same  path  as  that  of  the  1980-1981  ring. 
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Fig.  4.  Drifter  path,  inferred  center  path  (red),  and  orbit  ellipse 
(blue)  for  the  dnfter  in  the  1982-1983  ring.  Arrows  on  the  paths  are  at 
lO^lay  intervals,  while  the  ellipses  are  shown  at  15-day  intervals. 
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Fig.  5.  Swirl  (left  scale)  and  translation  velocities  (right  scale)  in¬ 
ferred  from  the  three  dnfters  in  the  1982-1983  ring.  The  east  and 
north  components  of  the  translation  velocity  are  depicted  in  red  and 
blue,  respectively. 

Figure  5  shows  the  translation  and  swirl  velocities  for  the 
1982-1983  ring.  The  swirls  are  slightly  less  than  50  cm/s, 
which  could  reflect  a  smaller  orbit  of  the  drifter.  The  dominant 
period  is  still  about  14  days.  Note  that  around  day  1,  1983  the 


swirl  velocities  decrease  significantly  and  show  some  evidence 
of  higher-frequency  components.  This  is  the  period  when  the 
drifter  moves  west  of-  94° W  and  executes  the  deformed, 
northwest-oriented  loop.  As  with  the  earlier  ring,  this  is  inter¬ 
preted  as  evidence  for,  topographic  aiid/or  slope  circulation 
interactions. 

Figure  2  indicates  that  west  of  94°W,  the  three  drifters  in 
the  1980-1981  ring  show  significant  divergence  in  the  calcula¬ 
tion  of  the  inferred  center.  Furthermore,  drifters  from  both 
rings  showed  anomalous  path  characteristics  in  this  region.  In 
particular,  1598  and  1600  appeared  to  become  entrained  in 
smaller,  anticyclonic  eddies,  while  1599  and  3374  exhibited 
warped  and  elongated  ellipses.  This  has  been  attributed  to 
interaction  between  the  ring  and  the  continental  slope  and/or 
slope  circulation  further  to  the  west. 

Figure  6  focuses  on  this  time  interval.  Figure  6a  shows  the 
paths  of  ail  four  drifters  during  the  period  in  question.  It  is 
clear  from  this  that  all  four  drifters  were  under  the  influence  of 
different  dynamical  processes  here  than  those  they  experienced 
to  the  east. 
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Fig.  66 

Fig.  6.  Paths  of  drifters  1598,  1599  and  1600  in  the  1980-1981 
ring  and  3374  in  the  1982-1983  ring  for  the  period  when  they  were 
interacting  with  the  slope  topography  and/or  circulation.  The  arrows 
are  at  10-day  intervals,  (a)  Drifter  paths.  (6)  Inferred  centers  of  flow. 

Figure  66  compares  the  movement  of  the  inferred  centers. 
In  interpreting  this  figure,  it  should  be  emphasized  that  the 
center  calculation  is  greatly  complicated  by  the  changes  tn 
curvature  of  the  paths.  The  paradigm  interprets  this  as  a 
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change  from  an  anticyclonic  to  cyclonic  structure  with  a 
consequent  change  in  the  location  of  the  center.  No  doubt  the 
real  dynamics  during  this  period  are  much  more  complicated 
than  the  rather  simple  scenario  available  from  the  kinematic 
model.  The  geographic  variation  in  the  paths  clearly  is  much 
larger  than  that  seen  east  of  94°W  or,  as  will  be  seen  shortly, 
than  that  found  along  the  slope.  This  suggests  that  each  of 
these  drifters  was  either  briefly  pulled  out  of  the  parent  ring  or 
its  orbits  deformed  through  interaction  with  the  continental 
slope  topography  or  circulation.  This  hypothesis  accounts  for 
the  deformed  orbits  and  the  brief  occurrence  of  high  fre¬ 
quencies  in  the  swirl  velocities. 

3.2.  Simulations 

Attention  is  now  turned  to  the  simulated  data.  There  are 
two  issues  involved  in  the  utilization  and  interpretation  of  this 
data.  First,  how  consistent  is  the  paradigm  in  recovering  ring 
kinematics  which  should  be  independent  of  the  orbit?  This 
issue  is  presently  being  addressed.  Addressed  here  is  the  issue 
of  establishing  how  well  the  GCM  agrees  with  the  observa¬ 
tions.  Speciiically,  the  center  path,  orbit  ellipses,  and  trans¬ 
lation  velocity  as  inferred  from  the  simulated  data  ate  com¬ 
pared  with  the  same  properties  as  determined  from  the  obser¬ 
vations. 
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Figure  7  shows  the  path  of  one  of  several  simulated  drifters 
from  the  HT  GCM  along  with  the  orbit  ellipses.  Figure  la 
depicts  the  results  using  the  ring  center  as  inferred  from  the 
drifter  orbit,  while  Figure  lb  repeats  the  ellipse  calculation 
using  the  position  of  the  maximum  displacement  anomaly  of 
the  interface.  The  choice  of  the  interface  displacement  as  rep¬ 
resentative  of  the  ring  center  is  arbitrary.  Presently,  we  are 
comparing  the  maximum  interface  displacement,  the  maxi¬ 
mum  surface  elevation  displacement  and  orbit-determined 
centers.  In  general,  all  are  within  30  km  of  each  other. 

The  most  significant  difference  between  the  orbit- 
determined  centers  and  interface  displacement  anomaly  occurs 
in  the  interval  between  9rW  and  93°W  where  the  orbit- 
inferred  ring  center  shows  a  significant  southern  loop.  This  is 
not  seen  in  the  path  of  the  maximum  interface  displacement 
anomaly.  At  93°W  the  orbit-inferred  center  executes  an  anti¬ 
cyclonic  loop  that,  again,  is  not  seen  in  the  displacement 
anomaly  path.  Elsewhere,  the  orbit-determined  path  is  consis¬ 
tently  S  to  25  km  south  of  the  displacement  anomaly  path. 

A  comparison  of  Figures  2, 4,  and  7  shows  that  the  simulat¬ 
ed  and  observed  centers  follow  the  same  general  path  through 
this  part  of  the  midgulf.  In  general,  the  simulated  center  path 
is  10  to  35  km  south  of  the  observed  path,  depending  upon 
which  of  the  two  simulated  center  paths  is  used.  There  is 
virtually  no  evidence  of  the  interaction  with  the  continental 
slope  topography  and  circulation  in  the  model  results  as  was 
seen  in  the  observations.  None  of  the  simulated  orbits  were 
deformed  or  appeared  to  be  pulled  out  temporarily  from  the 
ring.  Somewhat  surprisingly,  we  have  found  less  consistency  in 
the  ellipse  calculations  for  the  simulations  than  in  the  observa¬ 
tions.  A  number  of  the  ellipses  from  the  simulated  data  are 
quite  deformed.  They  also  show  frequent  reversals  of  the 
major  and  minor  axes.  Overall,  however,  as  the  ring  nears  the 
slope  region,  a  northeast  orientation  appears  to  develop  in  the 
orbit  ellipses. 

Figure  8  shows  the  translation  and  swirl  velocities  for  the 
simulation.  Figure  8a  gives  the  results  using  the  orbit  center, 
and  Figure  8h  uses  the  maximum  displacement  anomaly  as 
the  ring  center.  The  dominant  period  in  the  swirl  velocities  is 
about  24  days,  or  10  days  more  than  was  found  in  the  obser¬ 
vations.  Note  that  the  swirl  velocities,  as  calculated  from 
either  center,  are  of  comparable  magnitude,  about  30  cm/s. 


LONGITUDE  TIME  UULIAN  DAYS) 

Fig.  76  Fig.  80 

Fig.  7  Path  and  orbit  ellipses  (blue)  of  a  simulated  dnfter  from  Fig.  8.  Translation  and  swirl  velocities  as  inferred  from  the  HT 
the  HT  GCM.  Arrows  on  the  path  are  at  10-day  intervals,  while  GCM.  The  east  and  north  components  of  the  translation  velocity  are 

ellipses  are  shown  at  15-day  it  tervals.  Figure  7a  shows  the  ellipses  depicted  in  red  and  blue,  respectively.  The  scale  for  the  translation 

using  the  orbit-inferred  center  aj  the  flow  center,  and  Figure  76  uses  velocity  is  on  the  right.  Figure  8a  shows  results  using  the  orbit- 

the  maximum  interface  displacement  anomaly.  Both  centers  are  inferred  center,  while  Figure  86  {next  page)  shows  results  using  the 

shown  in  red.  maximum  displacement  anomaly  as  the  flow  center. 
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This  is  somewhat  less  than  the  observed  swirls,  but  this  could 
be  rectified  by  picking  a  simulated  drifter  in  a  greater  orbit. 
There  seems  to  be  more  difference  in  the  translation  velocities. 
Because  its  path  is  straighter,  the  displacement  anomaly  path 
has  less  variance.  The  mean  westward  component  in  both 
cases  of  about  5  cm/s  is  consistent  with  the  observations. 


LONGITUDE 

Fig.  9.  The  centers  for  the  1980-1981  ring,  the  1982-1983  ring, 
and  the  HT  GCM  simulation.  Orbit  1  (red)  refers  to  the  orbit- 
determined  center,  and  pressure  (blue)  refers  to  the  maximum  dis¬ 
placement  anomaly  path.  Arrows  depict  10-day  intervals. 

Figure  9  compares  the  inferred  center  movement  for  all  four 
drifters  along  with  the  inferred  center  of  the  simulated  dnfter 
and  the  displacement  anomaly  path  east  of  94°W.  This  is  the 
region  where  there  is  no  apparent  effect  of  slope  topography 
or  arculation  m  the  observations.  Of  course,  details  in  the 
paths  vary,  but  overall  the  agreement  is  quite  good.  The  maxi¬ 
mum  deviation  at  any  one  time  between  any  of  the  paths  is 


100  km,  which  is  only  two-thirds  the  diameter  of  a  typical 
ring.  For  most  of  the  paths,  the  orbit-determined  centers  are 
somewhat  south  of  the  displacement  anomaly  path.  However, 
as  the  paths  near  94°W,  the  observed  centers  move  north  of 
the  displacement  anomaly.  All  the  orbit-determined  center 
paths  show  considerably  more  structure  than  does  that  of  the 
displacement  anomaly. 

4.  Continental  Slope 

It  is  remarkable  that  both  observed  rings  and  the  simulated 
ring  impacted  the  continental  slope  in  the  same  area.  This  is 
seen  in  Figure  10,  which  shows  the  drifter  paths,  inferred 
center  paths,  and  orbit  ellipses  for  all  three  cases.  In  the  case 
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Fig.  10.  Drifter  paths,  inferred  centers  and  orbit  ellipses  for  (a)  the 
1981-1982  ring,  (6)  the  1982-1983  ring  (above);  (c)  the  simulated  ring, 
and  (4)  center  paths  along  the  continental  slope  (next  page).  Arrows 
on  the  paths  depict  lO-day  intervals,  while  the  ellipses  are  at  IS-day 
intervals. 

of  the  simulated  data  (Figure  10c),  the  orbit-inferred  center 
path  and  the  displacement  anomaly  path  are  virtually  identi¬ 
cal.  Consequently,  it  is  only  necessary  to  display  one  set  of 
calculations.  The  orbit-determined  centers  were  used  here. 

Figure  104  illustrates  how  close  the  movements  of  the  two 
observed  and  simulated  rings  are.  There  the  paths  of  the  in¬ 
ferred  centers  along  with  the  displacement  anomaly  path  have 
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been  superimposed.  Except  for  the  beginning  and  end  of  the 
paths,  the  maximum  deviation  is  70  km  with  a  rms  deviation 
of  12  km.  Moreover,  all  paths  show  the  same  general  charac* 
tefistics.  Both  the  observations  and  simulation  show  a  north¬ 
ward  migration  along  the  isobaths  but  with  some  eastward  or 
looping  component. 
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Fig.  10.  (com.)  Drifter  paths,  inferred  centers  and  orbit  ellipses  for 
(c)  the  simulated  ring,  and  (d)  center  paths  along  the  continental  slope. 
Arrows  on  the  paths  depict  10-day  intervals,  while  the  ellipses  are 
shown  at  15-day  intervals. 

The  northward  migration  along  the  slope  is  puzzling,  as 
conventional  theory  [Smith  and  O’Brien,  1983]  would  have 
topographic  beta  drive  the  rings  to  the  south.  Recently,  Smith 
[1986]  has  suggested  that  the  northward  migration  may  be 
due  to  boundary  effects.  An  alternate  explanation  was  offered 
by  Nakamoto  [1986],  who  found  northward  propagating  soli- 
ton  solutions  in  a  two-layer,  nonlinear,  quasi-geostrophic 
model. 

The  swirl  and  translation  velocity  time  series  are  shown  in 
Figure  11.  In  the  1980-1981  ring  the  dominant  period  is  15 
days,  almost  exactly  that  found  in  the  midgulf.  For  the  1982- 
1983  ring  the  dominant  period  is  about  11  days,  slightly  less 
than  that  found  in  the  midgulf.  The  simulated  ring  also  shows 


a  decrease  in  the  dominant-  period  to  about  15  days.  The 
magnitudes  of  the  observed  ring’s  swirl  velocities  are  about  40 
to  50  cm/s,  which  is  close  to  that  found  in  the  midgulf.  The 
swirl  velocity  for  the  simulated  ring  is  about  80  cm/s,  substan¬ 
tially  more  than  what  was  obtained  for  the  midgulf.  All  trans¬ 
lation  velocities  indicate  a  fairly  steady  northward  movement 
of  about  4  cm/s. 
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Fig.  11.  Swirl  (left  scale)  and  translation  velocities  (right  scale)  for 
the  (a)  1980-1981,  (b)  1982-1983,  and  (c)  simulated  rings  along  the 
continental  slope. 
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5.  Conclusions 

The  paths,  translation  and  swirl  velocities,  and  orbit  ellipses 
have  been  calculated  from  Lagrangian  observations  for  two 
large  anticyclonic.  Gulf  of  Mexico  rings  and  one  simulated 
ring.  Within  the  capabilities  of  the  data  base  and  the  analysis 
routine,  the  following  are  now  established. 

1.  The  paths  of  the  three  rings  across  the  central  Gulf  of 
Mexico  to  the  continental  slope  off  of  Mexico  are  virtually 
identical.  However,  not  all  Loop  Current  rings  will  follow  this 
path.  See  Lewis  and  Kirwan  [1987]. 

2.  The  intraring  variability  of  translation  and  swirl  veloci¬ 
ties  and  the  inferred  center  paths  for  the  1980-1981  ring  are 
essentially  the  same  as  the  inter-ring  variability  between  the 
1980-1981  and  1982-1983  rings. 

3.  For  the  observed  rings,  there  was  evidence  of  strong 
ellipticity  developing  as  these  rings  approached  the  conti¬ 
nental  slope.  The  axis  of  orientation  was  approximately 
NE-SW.  No  independent  data  on  ellipticity  were  available  for 
the  1982-1983  ring,  so  it  cannot  be  ruled  out  that  this  ring 
was  so  orientated.  Similar  but  less  pronounced  ellipticity  was 
seen  in  the  simulated  ring. 

4.  Both  observed  rings  showed  significant  distortion  of  the 
drifter  orbits  west  of  94°W.  This  is  attributed  to  interaction  of 
these  rings  with  the  continental  slope  topography  and/or  slope 
circulation  features.  Such  interaction  was  not  observed  in  the 
simulated  ring. 

5.  After  encountering  the  continental  slope  off  of  Mexico, 
both  observed  and  simulated  rings  retained  their  integrity  for 
at  least  several  months  while  they  migrated  to  the  north.  The 
track  of  the  center  of  both  was  very  close  to  the  1200-m 
isobath.  The  observed  swirl  speeds  were  slightly  changed  from 
the  values  obtained  from  the  midgulf.  The  simulated  ring  ex¬ 
hibited  a  substantial  increase  in  swirl  velocity.  For  both  the 
observed  and  simulated  data,  the  translation  speeds  were  com¬ 
parable  to  those  obtained  for  the  midgulf. 

There  are  a  couple  of  broader  implications  of  this  study 
pertinent  to  eddy-resolving  GCM’s.  Both  the  volume  of  data 
as  well  as  the  time  and  space  scale  resolution  of  simulated 
data  available  from  these  models  far  surpasses  the  typical  ob¬ 
servational  data  bases  used  for  verification.  Because  of  the 
different  time  and  space  scales  of  resolution  in  the  simulated 
and  observed  data  bases,  it  is  not  clear  how  to  establish  reli¬ 
able  criteria  for  statistical  comparisons.  This  study  suggests 
that  a  potent  alternative  is  to  compare  model  and  observed 
Lagrangian  kinematics.  Both  simulated  and  observed  La¬ 
grangian  data  sets  provide  comparable  space-time  resolution 
on  the  evolution  of  specific  circulation  features.  Moreover,  the 
assessment  of  the  prediction  of  such  features  is  a  more  strin¬ 
gent  test  of  a  model’s  predictive  capability  than  are  general 
statistics.  The  Lagrangian  data  also  provide  a  means  of  fine- 
tuning  model  parameters  such  as  layer  depth,  density  differ¬ 
ences,  etc.,  so  that  observed  and  simulated  swirl  velocities  can 
be  matched.  As  to  the  HT  model,  this  study  has  documented 
the  interesting  situation  that  the  prediction  improves  in  time, 
at  least  as  far  as  the  movement  of  the  ring  is  concerned. 

The  second  and  related  issue  is  the  use  of  Lagrangian  obser¬ 
vations  for  updating  prognostic  models.  In  character,  such 
data  are  similar  to  sea  surface  topographic  data  derived  from 
satellites  [see  Hurlburt,  1986;  Thompson,  1986;  Kindle,  1986]. 
The  Lagrangian  data  contain  information  on  the  motion  of 
representative  parcels  but  no  information  on  what  is  happen¬ 
ing  nearby  Satellite  topographic  data  return  information  just 


along  the  satellite  path  but  give  no  information  on  the  actual 
motion.  Supplethenting  the  latter  data  with  the  former  data  in 
eddy-resolving  general  circulation  models  could  significantly 
upgrade  their  predictive  capabilities. 

Appendix 

The  purpose  here  is  to  provide  some  details  on  the  inver¬ 
sion  of  (3)  and  (4)  to  obtain  the  time  series  for  the  elements  of 
Z^.  The  basis  for  this  inversion  is  a  Taylor  expansion  in  time 
about  the  instant  t*  of  the  velocity  vector.  From  the  left-hand 
side  of  (3)  and  (4)  one  obtains 

u(t)  =  liftj)  +  u'(fjXt  -  t*)  +  -  tt)V2 

+  «'"(t»Xt  -  MV6  (Al) 

tXt)  =  tXti)  -  '  .tjXt  -  ft)  +  -  tk)^/2 

+  t>"'(t»Xt  -  h)V6  (A2) 

for  the  interval  The  primes  represent  time  de¬ 

rivatives.  Now  each  of  the  terms  on  the  right-hand  side  of  (Al) 
and  (A2)  can  be  estimated  from  the  velocity  record.  For  exam¬ 
ple  u(tt)  and  i^tjt)  are  merely  the  velocities  at  time  The 
derivatives  can  be  estimated  by  a  variety  of  techniques;  here 
we  have  employed  centered  finite  differences. 

Now  the  analytic  solution,  (S)  and  (6),  can  be  expanded  in  a 
Taylor  series  as  well.  When  this  is  done  and  coefTicients  of  the 
appropriate  powers  in  this  expansion  and  (Al)  and  (A2)  are 
equated,  a  system  of  simultaneous  nonlinear  equations  is  ob¬ 
tained  for  each  time  interval.  With  the  subscript  k  suppressed, 
these  are 


X{a  +  d}+  Y{b  -  c)  -f  2Ur  *  2« 

(A3) 

X(b  +  c)  -f  Y{d  -«)•+•  IVt  =  2v 

(A4) 

Xlid  +  o)*  +  A’  -  c*]  +  2Yd(b  -  c)  =  4u' 

(A5) 

2Xd{b  -1-  c)  +  y[(d  -  fl)*  +  b^-  c*]  =.  4i;' 

(A6) 

X[{d  +  0)’  +  (A*  -  c*X3d  +  a)] 

+  Y(b  -  cX3d*  +  a*  +  A*  -  c*)  =  8u" 

(A7) 

X(b  +  cX3d^  +  o’  +  A^  -  e’) 

+  y[(d  -  fl)’  +  (A’  -  c’X3d  -  0)]  =  8v" 

(A8) 

A’Ko*  +  6*  -  c*  -t-  d*)*  +  4d*(a’  +  —  c*) 

+  4ad(a^  +  b^  —  c^  +  d*)] 

+  4Yd{b  -  cXo’  +b^-c^  +  d*)]  =  16u"'  (A9) 

4Xd{b  +  cXo*  +  6*  -  c*  +  d*)  +  y[(o*  +  b^-c^  +  d»)* 

+4d*(fl’+6*— c’)— 4<id(a’  +  6’ — c* + d*)]  =  16i;'"  (A  10) 

These  eight  equations  are  inverted  at  each  time  step  for  X,  Y, 
I/r>  ^bd  d.  To  our  surprise  it  seems  that  this  can  be 

done  analytically.  The  key  to  this  is  the  observation  that  the 
geometric  invariants  of  the  matrix  M^,  i.e.,  Tr  (Af)  and  det  (M), 
can  be  calculated  from  observations  without  knowing  a,  b,  or 
c,  a  priori.  After  3  years  of  inspection  it  was  seen  that 

det  (Af)  =  2(u"v"'  -  u-'VWv"  -  v‘u“)  =  (All) 

Tr  (Af)  =  (u'v'"  -  t}'u'"Wv"  -  v'u")  =  d  (A12) 

Insertion  of  Af^  and  d,  calculated  from  (All)  and  (A12), 
respectively,  into  (A5>-<A10)  significantly  simolifies  the  latter 
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expressions.  Considerable  routine  algebra  then  yields 


X  =  8[u'(M^  -F  2d^)  ^ 

(A13) 

Y  =  8[i''(M^  +  2d^)  -  2v"‘2/M* 

(A14) 

a=  +  Ki  +  Ki(dl2)}/N 

(A15) 

b  =  -1-  Hj  -F  Hi{d/2)}/N 

(A16) 

c  =  -{G,M^  -  Gj  -F  Gi{dJ2)}/N 

(AI7) 

where 

N  =  8(«V"  - 

G,  =  2{u'^  +  v'^); 

Gi  =  8(u'^  +  v"^)-, 

Gj  =  16(u'u''  +  v'v"); 

H,  =  2(u'^  -  v‘^y. 

Hi  =  8(u"*  -  y"^); 

Hi  —  I6(u'u''  —  u"o'); 

Ki  =  2m'u'; 

Ki  =  8u"v"; 

K,  =  8(u'y”  +  u"v'). 

Incidentally,  the  parameter  N  is  directly  related  to  the  angu¬ 
lar  momentum  per  unit  mass  of  the  orbit  L: 

L  =  -4N/M*  (A  18) 

The  final  step  in  the  calculations  is  to  determine  the  trans¬ 
lation  velocities  from  (A3)  and  (A4).  This  is  trivial,  since  all 
other  terms  are  now  known. 
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Notes  on  the  Cluster  Method  for  Interpreting  Relative  Motions 
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IvnooucnoN 

Esismates  of  the  horizontal  spatial  derivath'cs  of  the  near- 
surface,  horizontal  velocity  are  of  obvious  importance  in  a 
wide  range  of  scales  in  ocean  dynamics.  They  are  notoriously 
didicult  to  obtain.  Most  estimates  have  come  from  rdative 
motion  analysis  of  drifter  dusters.  For  an  older  summary  of 
results,  see  Kirwen  £1975}.  Fahrbach  et  oL  [199^  hav-e  recent¬ 
ly  applied  this  technique  to  drifter  data  from  the  equatorial 
Atlantic. 

Sanderson  et  at.  [1988]  have  pointed  out  a  statistical  bias  in 
using  this  technique  if  the  number  of  drifters  in  the  duster  is 
small.  Of  course,  the  statistical  significance  of  modd  parame¬ 
ter  estimates,  as  well  as  the  effect  of  measurement  errors  in 
estimating  the  velodty  gradient  tensor,  is  important  There  ». 
however,  another  problem  with  this  technique.  This  is  the  tadt 
assumption  that  the  center  of  mass  of  the  cluster  corresponds 
to  the  center  of  the  flow  field  about  which  a  Taylor's  ex¬ 
pansion  of  the  velodty  field  can  be  made. 

The  purpose  here  is  twofold.  The  first  is  to  show  that  with 
the  cluster  paradigm,  a  deviation  of  the  cluster  center  from  a 
flow  center  can  result  in  serious  biases  in  gradient  estimates. 
The  other  is  to  propose  an  alternative  to  the  cluster  model 
that  avoids  the  center  bias  problem.  The  new  approach  pro¬ 
vides  independent  estimates  of  the  velodty  gradient  estimates; 
however,  it  requires  solving  a  nonlinear  inverse  problem. 

Theory 

The  first  part  of  the  cluster  paradigm  is  given  by 

('■/  -  =  “/  =  U/  +  I  r./  -H  r/7A  (1) 

* 

Here  A  is  the  time  difference  between  fixes,  is  the  a  compo¬ 
nent  of  the  position  vector  from  the  cluster  center  of  mass  to 
the  drifter  path,  is  its  original  or  previous  position  vector, 
is  the  velodty  gradient,  and  U/  is  a  steady  mean  flow. 
Figure  la  is  a  cartoon  of  the  general  geometry.  As  has  been 
noted  by  many  previous  investigators,  (1)  can  be  obtained  by  a 
Taylor's  expansion  of  the  velodty  field. 

The  first  term  on  the  right-hand  side  of  (1)  is  a  uniform 
translation  of  the  center,  and  the  second  term  is  a  local  swirl 
induced  by  a  homogeneous  deformation  field.  The  last  term  in 
(1)  is  a  residual,  the  minimization  of  which  yields  estimates  of 
and  For  the  latter  operation  it  is  necessary  to  intro- 
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duoe  the  duster  average  of  property  h' 

<*^  =  £A'/A7  12) 

where  \*  is  the  statistical  degrees  of  freedom  in  the  duster 
average  and  the  sum  is  over  all  drifters  in  the  duster  [see 
Sanderson  et  aL.  1988].  Note  that  by  definition 

<0=0  (3) 

This  is  the  second  part  of  the  duster  paradigm.  It  is  stressed 
that  (3)  significantly  constrains  the  basic  assumption  and  va¬ 
lidity  of  a  two-term  Taylor  expansion. 

By  standard  minimization  augmented  with  (3).  the  least 
squares  estimates  of  V/  and  are 


1// =<«.'> 

(4) 

(5) 

7 

(6) 

is  the  inverse  of  the  statistical  moment  of  inertia  matrix  of  the 
cluster.  Details  and  variations  of  (4)  and  (S)  are  given  by  Molin- 
art  and  Kirwan  [1975],  Okubo  and  Ebbesmeyer  [1976],  and 
Okubo  et  al.il976\. 

This  model  is  critically  dependent  upon  (1)  a  translation 
velodty  U/  common  to  all  drifters,  (2)  a  velodty  gradient 
tensor  common  to  all  drifters,  Le.,  a  homogeneous  defor¬ 
mation  field,  and  (3)  equation  (3). 

In  order  to  assess  the  role  these  assumptions  play,  consider 
a  more  general  kinematic  flow  model: 

“/  =^/  +  lL  +  V'/A  (7) 

b 

Here  is  the  position  vector  to  the  drifter  path  from  a 
physically  relevant  local  origin.  This  is  called  the  flow  center 
to  distinguish  it  from  the  cluter  center.  See  Figure  la  for  the 
geometric  relation  between  x/  and  r/.  The  latter  depends 
upon  the  geometric  characteristics  of  the  drifter  array;  the 
former  depends  on  true  flow  field  characteristics.  This  model 
differs  from  (I)  in  that  it  allows  for  a  separate  Taylor's  ex¬ 
pansion  for  each  drifter.  Thus  there  are  distinct  translation 
velocites  V/,  as  well  as  distinct  velodty  gradients  for 
each  drifter.  Also  in  contrast  to  (3) 

<x/>?4  0  (8) 

It  should  be  clear,  however,  that  if  the  three  assumptions  listed 
in  the  preceding  paragraph  are  imposed  on  (7),  it  will  reduce 
tod). 
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vdodiy  and  is  of  the  order  of  the  true  velocity  gradient  times 
tbe  distanoe  betiween  tbe.two  centers.  For  small  area  clusters, 
this  dispiaoement  could  «dl  be  of  the  duster  horizontal  scale, 
and  so  the  duster  determined  translation  ydodty  could  be  of 
tbe  order  of  the  swirl  velodty.  In  tbe  case  of  tbe  vdodty 
gradient  estimate,  there  is  no  a-prion  reason  to  expect 
to  be  small:  hence  tbe  “borrectioh**  term  in  (13)  could 
be  of  tim  order  of  the  correct  value.  As  an  aside.  (13)  is  recog- 
mzed  as  an  extension  of  Pappas'  theorem  from  mechanics. 

Case  //.  Tbe  translation  vdodty  varies  from  drifter  to 
drifter:  however,  assumptions  2  and  3  hold  (see  ITgure  Ic). 
This  case  would  indude  tbe  situation  where  tbe  drifters  are 
deployed  in  a  jet  such  v  tbe  Gulf  Stream.  Perhaps,  fortu¬ 
itously.  tbe  duster  centd  initially  coincides  with  the  flow 
center,  say.  tbe  axis  of  tbe  jet.  Tbe  duster  parameter  estimates 
now  are 


Fig.  I.  Cartoon  showing  tbe  essential  geometry,  {a)  Basic  vectocs 
z.  r.  and  x  for  tbe  present  configuration,  (b)  Displacements  for  case  I. 
This  could  result  from  equal  shear  and  nonnal  deformation  and  anti- 
cjdonic  roution.  The  solid  line  indicates  the  present  configuration, 
and  tbe  dashed  line  indicates  the  previous  configuration,  (c)  Displace- 
ments  for  case  II.  Initially,  tbe  cluster  and  flow  centers  coincide,  but 
in  the  present  (solid  line)  configuration  the  flow  center  has  moved 
away  from  the  duster  center.  All  panels  depict  tbe  movement  and 
deformation  of  tbe  material  lines  connecting  the  drifters  as  well  as  a 
hipotbetical  array  initially  symmetric  about  tbe  flow  center. 

Given  only  cluster  data,  the  observables  are  V/  and  r/. 
We  now  examine  the  effect  of  naively  applying  the  cluster 
paradigm  to  this  more  general  situation.  Taking  the  cluster 
center  as  the  local  origin  and  letting 

r/  =  r/-x/  (9) 

be  the  displacement  vector  from  the  cluster  center  to  flow 
center,  it  is  deduced  from  (4),  (5),  (7),  and  (9)  that 

f//  =  <F/>-Hl<G^'(r/-z.')>  (10) 

» 

=  I  «»,«»'/'•,>  +  I  W>  (1 1) 

T  a 

These  last  two  equations  show  that  in  general,  the  cluster 
model  parameters  are  weighted  averages  of  the  flow  model.  In 
order  for  (10)  and  (11)  to  reduce  to  (4)  and  (5),  all  three  as¬ 
sumptions  given  above  must  be  made. 

The  following  special  cases  illustrate  the  tenuous  character 
of  *he  cluster  model. 

Case  1.  There  is  no  translation  (VJ  s  0),  and  the  defor¬ 
mation  field  is  homogeneous;  i.e.,  is  the  same  for  all 
drifters  (G^'  =  G^).  An  example  of  this  situation  would  be 
when  the  drifters  are  embedded  in  a  flow  field  characterized  by 
no  translation  and  homogeneous  vorticity  and  deformation 
fields  (see  Figure  lb).  It  is  seen  immediately  from  (10)  and  (1 1) 
that  the  cluster  model  will  yield  bogus  parameter  estimates  of 

I//  =  -  I  (12) 

b 

=  (13) 

r  ^ 

The  bogus  translation  velocity  has  the  structure  of  a  swirl 


Vs'  =  <K'>  (14) 

H-  =  +  I  (15) 

r 

If  the  flow  center  is  indeed  tbe  jet  axis,  then  the  cluster 
model  will  yidd  a  translation  velocity  for  the  center  somewhat 
less  than  the  maximum  velocity  of  the  jet  as  indicated  by  (14). 
Moreover,  if  tbe  jet  velocity  is  symmetric  about  the  flow  axis, 
then  tbe  cluster  gradient  estimates  of  the  cluster  scale  velocity 
gradients  will  include  an  alias  from  the  larger-scale  shear 
across  tbe  jet  as  exemplified  in  tbe  inhomogeneous  nature  of 
V/  (the  second  right-hand  side  term  in  (13)).  In  general,  the 
cluster  model  may  not  be  able  to  distinguish  between  small- 
scale  and  large-scale  shears. 

These  two  cases  highlight  a  basic  diiliculty  of  the  cluster 
method.  If  the  area  covered  by  the  cluster  is  small,  there  is  a 
high  likelihood  that  there  will  be  a  significant  displacement 
between  tbe  cluster  and  flow  centers.  This  is  where  case  I 
applies.  This  can  be  alleviated  by  increasing  tbe  area  extent  of 
the  cluster  (with  or  without  additional  drifters),  but  then  the 
large-scale  shear  may  defeat  the  cluster  model,  as  is  indicated 
by  the  analysis  in  case  II. 

Case  in.  Only  assumption  2  holds.  In  this  case,  the  clus¬ 
ter  parameter  estimates  will  superpose  (12)  and  (14)  for  U/ 
and  (13)  and  (15)  for 

It  is  emphasized  that  the  problems  indicated  above  cannot 
be  rectified  by  increasing  the  number  of  drifters  and/or  the 
cluster  area  and/or  the  position  fix  accuracy  and  frequency. 
The  basic  issue  here  is  the  weakness  of  the  physics  implied  by 
the  approach,  not  the  accuracy  or  statistical  significance  of  the 
measurements. 

Discussion 

The  preceding  analysis  suggests  that  the  cluster  method  for 
estimating  velocity  gradients  should  be  used  with  extreme  cau¬ 
tion.  Other  observations  should  be  used  to  establish  its  appli¬ 
cability  in  specific  experiments. 

There  is  an  alternative  to  the  cluster  paradigm  for  extract¬ 
ing  velocity  gradient  information  from  path  data.  The  alter¬ 
native  makes  use  of  a  kinematic  model  first  proposed  by 
Okubo  [1970]  and  recently  discussed  by  Perry  and  Chong 
[1987].  This  approach  makes  use  of  the  fact  that  (7)  (and 
hence  (1))  can  be  readily  integrated  to  obtain  particle  paths.  If 
the  deformation  field  is  homogeneous  and  stationary,  then  a 
single  particle  path  contains  the  same  information  on  the  de¬ 
formation  as  does  a  hypothetical  cluster  of  nearby  particles. 
The  only  characteristic  difference  in  the  paths  is  their  starting 
positions. 

To  illustrate  how  the  deformation  information  is  contained 
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in  the  path  data,  it  is  necessary  to  calculate  the  path.  Neglect¬ 
ing  the  last  term  in  (7)  and  assuming  that  V/  is  steady  in  time, 
standard  techniques  yield 

H-.W  -  -  O  +  (16) 

-  W'i'-CO  =  -  O  +  xAO  (17) 

Here  is  the  absolute  position  at  time  r,  and  H^Hto)  is 
the  initial  position  of  the  flow  origin  at  time  i,  of  drifter  p. 

The  other  quantities  in  ( 16)  and  (17)  are 

x/{t)  =  {exp  [m.'ft  -  OJLXj'd'  +  O')  + 

-{cxp[m.'(t-tMXAG'-in+X,>’G,^']}/2I’’ 

(18) 

x/(f)  =  {exp  Cm.'(f  -  gjCJf/  Gj/  +  -  0^]} 

+  {exp[in_'(t-0][-A-/G,/+Jfi'(/'+G^]}/2/' 

(19) 


m,,.'  =  :TriG^n±n/2  (20) 

G'  =  G,/-G„'  (21) 

/'  =  {(Gn^  +  (G,/  +  G^-^"  -  (G,/  -  Gj.')"}‘«  (22) 

2f/  =  X/do)  (23) 


Given  the  path  data,  i.e.,  the  left-hand  side  of  (16)  and  (17) 
for  a  few  intervals  of  time,  the  velocity  gradient  G,/  and  trans¬ 
lation  velocity  V/  can  be  estimated.  This,  however,  is  a  prob¬ 
lem  in  nonlinear  parameter  estimation  and  hence  of  more 
fundamental  mathematical  interest  than  purely  linear  statis¬ 
tical  models.  Kinvan  et  al.  [1988]  have  provided  one  method 
of  obtaining  these  estimates.  Presumably,  more  efficient  meth¬ 
ods  are  available. 

It  is  stressed  that  the  estimates  of  the  velocity  gradient  and 
translation  velocity  apply  only  to  drifter  p.  Other  drifters,  if 
present,  would  supply  independent  estimates.  This  would 
allow  statistical  assessments  of  the  uniformity  of  the  flow  field 
parameters. 

A  cautionary  note  on  application  of  the  proposed  technique 
should  be  made.  First,  gradient  estimates  from  this  approach 


are  truly  Lagrangian,  and  m  there  may  be  some  difficulty 
relating  them  to  Eulerian  estimates.  Second,  there  are  certain 
pathological,  situations  where  the  deformation  field  will  pror 
duce  variability  in  only  one  of  the  velocity  components.  In 
these  ciues  it  may  not  be  possible  to  terover  all  the  gradient 
information.  See  Okubo  [1970]  and  Perry  and  Chong  [1987] 
for  a  complete  discussion  of  patholo^cal  situations  in  critical 
points  in  flow  fields. 
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